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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS: 

This  report  for  DAAG29-83-K-0061  supplements  K-0140  which  covered  the  first  twenty- 
one  months  of  ARO  support  for  our  Studies  of  Silicon-Refractory  Metal  Interfaces.  The 
emphasis  here  is  on  the  final  15  months  of  that  program.  It  was  funded  under  two  distinct 
grants  because  of  the  transfer  of  the  PI  from  the  University  of  Wisconsin  to  the  Univer¬ 
sity  of  Minnesota. 

overall  goal  in  this  project  was  to  examine  the  electronic  and  morphological 
interactions  which  occur  at  interfaces  involving  Si  and  metal  overlayers.  During  the 
course  of  this  project,  we  examined  overlayers  of  Cr,  Ti,  Sm,  Au,  Ca,  Ce,  and  Cu  using 
photoelectron  spectroscopy.  Auger  spectroscopy  and  LEED.  A  new  technique  involving 
angle-resolved  Auger  profiling  was  developed  to  determine  the  atomic  density  profile 
for  interfaces,  as  discussed  in  the  enclosed  papers.  >Its  utility  was  demonstrated 
through  studies  of  Cu/Si,  Au/Si,  and  Ce/Si.  Exhaustive  use  of  synchrotron  radiation 
photoemission  has  lead  to  a  detailed  understanding  of  the  electronic  bonding  of  silicides 
and  metal/silicon  interfaces.  Collaborative  studies  with  theorists  gave  rise  to  modeling 
of  interfaces  and  calculation  of  electronic  energy  states  for  ordered  silicides.  Gf-'i  •  h  , 
\  <  V  v  1-.  ■»  K\di,  G 

MAJOR  BREAKTHROUGHS  INCLUDE: 

■  A  systematic  evaluation  of  bonding  .in  silicides,  showing  how  the  metal -d/si 1 icon-p 
bonds  will  change  as  a  function  of  transition  metal  variation.  This  work  built  on  the 
work  by  others  for  the  Pd  and  Ni  silicides  and  created  a  broader-based  picture.  Experi¬ 
mental  studies  of  bulk  silicides  and  theoretical  and  band  structure  calculations  made 
this  panoramic  view  possible. 

Characterization  of  the  onset  of  reaction  on  silicon  surfaces,  showing  that  the 
details  vary  significantly  with  the  overlayer  character  but,  at  the  same  time,  that  all 
of  the  transition  metals  induce  reactive  disruption  of  the  surface.  For  Cr/Si ,  we  showed 
that  the  onset  of  reaction  was  delayed  until  a  critical  coverage  was  reached,  corres¬ 
ponding  to  patch  formation  and  subsequent  disruption.  For  Sm/Si  and  Ca/Si ,  the  reaction 
was  triggered  at  the  very  lowest  coverages.  For  Ce/Si,  an  even  more  complicated  mor¬ 
phology  was  observed  and  modeled.  At  low  coverage  it  was  found  that  clusters  formed 
but  that  for  coverages  approaching  0.6  monolayers  the  clusters  became  unstable  and  pro¬ 
vided  the  energy  necessary  to  disrupt  the  surface.  Such  an  energy  source  has  not  been 
demonstrated  for  Si -based  surfaces.  The  resulting  interface  was  then  highly  hetero¬ 
geneous.  To  model  it,  we  refined  a  method  of  determining  the  coverage  dependence  of 
the  three  Si  species  which  were  observed  to  form,  corresponding  to  substrate,  reacted, 
and  floating  Si. 

Extensive  modeling  of  interface  development  has  been  undertaken,  using  high  resolution 
photoemission.  With  this  technique,  we  are  now  able  to  fingerprint  the  species  which  form 
and  identify  the  chemical  profile  of  the  interface.  This  work  is  then  supplemented  by  the 
angle-resolved  Auger  studies  which  give  the  atomic  profile  of  the  interface.  Inportant 
information  regarding  in-  and  out-diffusion  of  metal  and  Si,  the  width  of  the  interface, 
the  dynamics  of  the  interface,  and  the  reactivity  of  the  interface  can  then  be  obtained. 
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MAJOR  BREAKTHROUGHS  (oont ’d) : 

The  work  accomplished  under  the  first  three  years  of  the  ARO  program  set  the  stage 
for  a  second  program  which  will  seek  to  understand  more  complex  systems,  including 
reactions  of  metal  overlayers  on  silicides,  silicon  overlayers  on  metal  surfaces,  and 
the  effect  of  thin  native  oxides. 
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Interface  catalytic  effect:  Cr  at  the  Si(l  1 1)-Au  interface 

A.  Franciosi  and  i.  H.  Weaver 

Department  »i  Chemical  Engineering  and  Materials  Science,  University  of  Minnesota. 

Minneapolis.  Minnesota  55455 

D.  G.  O'Neill 

Synchrotron  Radiation  Center,  University  of  Wisconsin — Madison. 

Stoughton.  Wisconsin  53589 
(Received  19  July  1983) 

Synchrotron  radiation  phoioemission  studies  of  the  effect  of  Cr  interlayers  on  Si(  1 1 1  )-Au  interface  reac¬ 
tion  show  that  Cr  concentrations  below  IxiO1’  atotns/cm-  retard  Si-Au  intermixing,  concentrations 
between  1  and  7  5x10'*  aioms/cm!  promote  Si-Au  intermixing,  and  concentrations  in  excess  of  8*I015 
j(oms/cm:  sharply  reduce  intermixing.  These  variations  are  shown  to  depend  on  the  three  formation 
stages  of  the  Si-Cr  junction  Cr  itself  is  shown  only  to  be  indirectly  involved  in  the  Si-Au  reaction  and  Si  is 
to  be  the  only  moving  species. 


The  driving  forces  responsible  for  atomic  interdiffusion  at  analyzer!  was  0.3-0. 5  eV  Cr  and  Au  were  evaporated  in 

Si-metal  interfaces  have  been  the  subject  of  intense  discus-  situ  from  W  coils  onto  cleaved  n-type  Si  crystals  (P-doped, 

sion.1"1  To  understand  this  atomic  interdiffusion.  we  must  !015  cm-1)  and  the  overlayer  thickness  was  measured  by  a 

identify  the  parameters  that  control  interdiffusion  kinetics  quartz-crystal  monitor  ( (»c,  ”  I  A  —  1.1  monolayers; 

and  the  junction  profile,  and  we  must  determine  the  rela-  9Xu-l  A -0  8  monolayers:  I  monolayer  -7.6xl014 

ttonshtps  that  exist  between  these  parameters  and  the  chem-  atoms/cm;)  Cr  depositions  of  0.5-15  A  were  made  im- 

ical  activity  of  the  species  involved  in  the  interface  forma-  mediately  before  a  series  of  Au  depositions  (maximum  Au 

tion  process.  ''1  coverage  of  50-70  A  for  each  Cr  interlayer  thickness!.  All 

In  order  to  better  understand  interdiffusion,  we  per-  experiments  were  performed  at  pressures  of  5  x  I0~"  Ton- 

formed  an  investigation  of  the  effect  of  Cr  interlayers  on  <  «6  x  10' 10  Torr  during  evaporation!  and  at  room  tempera- 

the  atomic  interdiffusion  of  Si  and  Au  at  the  Sill  ID-  ture. 

(2  x  U-Au  interface,  with  special  emphasis  on  the  chemical  In  Fig.  1  we  show  representative  photoelectron  energy 

aspects  of  interdiffusion.  distribution  curves  (EDC’s)  for  the  valence  band  of  the  Si- 

In  this  paper,  we  show  that  for  Cr  coverages  from  2  to  9  Cr-Au  interface  at  a  fixed  Au  coverage  of  20  A  as  a  func- 

A  the  interlayer  promotes  the  Si-Au  intermixing,  while  the  tion  of  Cr  interlayer  thickness.14  The  topmost  spectrum  cor-, 

interdiffusion  is  reduced  for  Cr  coverages  above  9  A—H0  responds  to  Si(lll)-Au  without  Cr,  and  the  bottom-mtist 

monolayers!  or  for  coverages  less  than  one  monolayer.  shows  the  valence  bands  for  a  20-A  Au  film  on  an  inert 

This  strikingly  nonmonotonic  behavior  rules  out  simple  substrate. 15  The  differences  between  these  two.  which  re- 

“diffusion  barrier”  effects  and  is  related,  instead,  to  the  veal  the  formation  of  the  Si-Au  interface,  have  been  ex- 

three  stages  of  Si-Cr  reaction.  Each  stage  corresponds  to  a  plained  as  modifications  of  Au5d-derived  density  of  states 

different  morphology  or  microscopic  arrangement  of  the  Si  features  caused  by  Si  atoms  in  the  Au  matrix,  mainly 

atoms  at  the  surface  and.  hence,  a  different  energy  content  through  the  change  of  the  d-d  overlap.1214 16 17  Taking 

of  the-Si-Si  bonds.  such  differences17  as  a  fingerprint  of  Si-Au  reaction,  the  sys- 

Our  results  show  that  the  interlayer  Cr  atoms  do  not  ap-  tematics  in  Fig.  1  show  that  a  Si-Au  reaction  can  occur  in 

pear  directly  involved  int  the  Si-Au  reaction.  Instead,  since  the  presence  of  a  Cr  interlayer,  but  that  the  intermixing  is 

the  Cr  interlayer  affects  the  formation  of  the  Si-Au  junction  sharply  reduced  when  the  interlayer  thickness  increases 

and  its  concentration  modulates  the  Si-Au  reaction,  it  seems  above  9  A. 

natural  to  define  the  role  of  the  Cr  atoms  as  that  of  an  “in-  The  presence  of  Si  atoms  in  the  Au  matrix  is  clearly  indi- 

terface  catalyst”  for  the  Si-Au  reaction.  cated  by  the  topmost  spectra  of  Fig.  1,  but  the  valence-band 

The  Si-Cr-Au  system  was  chosen  because  the  Si-Cr  and  spectra  show  no  evidence  of  Cr  at  the  surface.  To  deter- 

the  Si-Au  interfaces  are  well  characterized l0"u  and  neither  mine  whether  Cr  intermixes  with  the  Au  overlayer,  we 

gives  rise  to  island  formation  at  room  temperature.  The  measured  the  integrated  Cr  derived  3p  core  emission  as  a 

relative  surface  concentration  of  the  different  species  and  function  of  Au  coverage  for  two  different  Cr  interlayer 

the  evolution  of  the  chemical  bonding  was  followed  sys-  thicknesses,  as  shown  at  the  top  of  Fig.  2.  These  results 

tematically  through  synchrotron  photoemission  from  show  an  exponential  attenuation  of  the  Cr  emission.  Corn- 

valence  and  core  electronic  states,  both  as  a  function  of  Au  parison  with  a  simple  exponential  attenuation  calculated 

coverage  and  Cr  interlayer  thickness.  The  experiments  were  with  an  escape  depth  of  5  A  (dashed  line)  indicates  Si  out¬ 
performed  at  the  University  of  Wisconsin  Synchrotron  Radi-  diffusion  through  the  interlayer  into  the  Au  layer.  This  Si 

ation  Center  with  a  “Grasshopper”  grazing-incidence  mono-  outdiffusion  is  quantitatively  more  important  for  an  inter- 

chromator  for  40  «  h  v  *5 140  eV.  The  photoelectrons  were  layer  thickness  of  2  A  than  for  an  interlayer  thickness  of  8 

analyzed  by  a  double-pass  cylindrical  mirror  energy  A.  There  is  no  evidence  of  Cr  outdilTusion  into  the  Au 

analyzer,  and  the  overall  resolution  (monochromator  +  film. 
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FIG.  I  Valence-band  EDC's  of  the  Si-Cr-Au  system  as  a  func¬ 
tion  of  the  Cr  interlayer  thickness  at  a  fixed  Au  coverage  of  20  A. 
The  variation  of  the  energy  separation  of  the  Au  id  mam  features 
•  vertical  linesl  emphasizes  the  Sl-Au  intermixing  takes  place  in  the 
presence  of  the  interlayer  for  Cr  coveraie  below  9  A  but  is  sharply 
reduced  for  higher  interlayer  thickness 


intensity  of  the  Crip  cores  as  a  function  of  Au  coverage  on  the  Si- 
Cr  surface.  We  show  experimental  data  lor  Cr  interlayer 
thicknesses  of  2  and  8  A.  The  dashed  line  represents  the  theoretical 
result  for  exponential  attenuation  of  the  ip  based  on  an  electron  es¬ 
cape  depth  i  *5  A  Midsecuon  binding  energy  of  (he  Crip  cores 
as  a  function  of  Cr  coverage  for  Cr-Si<  1 1 11-12  *  1 1  The  Si-Cr  reac¬ 
tion  corresponds  io  a  chemical  shift  of  —0  3  eV  for  the  Cr3pcores 
Bottom  variation  of  the  ip  core  binding  energy  for  the  Cr  atoms  in 
the  interlayer  as  a  function  of  Au  coverage  Within  experimental 
uncertainty  ihe  average  chemical  environment  of  the  Cr  atoms 
remains  the  same  alter  ihe  Si-Au  mterdiffusion 


In  a  series  of  previous  papers  we  showed  that  the  deposi¬ 
tion  of  Cr  Into  the  Sill  11)  surface  produces  a  Si-Cr  mixed 
surface  phase  of  variable  composition."  When  the  Si  .r- 
Au  junction  is  formed,  we  now  see  that  St  diffuses  into  the 
Au  matrix  and.  in  principle,  this  might  imply  Si  depletion  in 
the  Si-Cr  region  However,  the  results  shown  in  Fig.  2  indi¬ 
cate  that  this  is  not  the  case.  The  Si-Cr  reaction  which  oc¬ 
curs  during  deposition  of  Cr  on  Si<  1 1 1 )  produces  a  chemical 
shift  as  large  as  0.3  eV  for  the  Cr3p  cores  •  midsection  of 
Fig.  2).  but  there  is  no  variation  of  this  cnemtcal  shift  when 
Au  is  deposited  onto  the  surface  •  lowest  section!  Hence, 
within  experimental  uncertainty  of  —01  eV.  the  average 
local  chemical  environment  of  Cr  atoms  in  the  Si-Cr  phase 
remains  the  same  after  the  Si-Au  mterdiffusion.  This  im¬ 
plies  that  as  Si  enters  the  Au  overlayer  a  corresponding 


amount  must  enter  the  Si-Cr  phase  to  maintain  the  local 
stoichiometry  The  net  results  of  this  diffusion  is  transport 
of  Si  atoms  through  a  stable  Si-Cr  interface  phase. 

In  Fig.  3  we  show  very  important  results  which  .'rmon- 
strate  that  a  picture  in  which  ihe  Si-Cr  phase  acts  as  a  dif¬ 
fusion  barrier  lor  atomic  mterdiffusion  of  Si  and  Au  is 
inadequate  Bv  measuring  the  attenuation  jf  the  Si  Ip  core 
emission  as  a  : -.action  of  Au  overlayer  tor  different  Si-Cr 
interlayer  thick  "esses,  'e  .an  sr.ow  tnat  the  attenuation  is 
not  proportion.-.'  to  :  —  nicicness  of  me  Cr  barrier.''  In 
Fig.  3  the  enis.-m  .rttensitv  .;.  a  giver.  Au  coverage  and  in¬ 
terlayer  thsckr.e  -  -  -  — ia.  e.vj  'o  ;he  initial  Si  2r  emission 
from  the  n-  onu-.  p  •*  c.  meanson  results  lor  the 
Si 1 1 11 1  -  A  u  interface  uittrou;  Cr  deposition  are  ...»o  shown 
(dot-dashed  line  from  Ref.  Ui  The  best  representation  of 
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FIG  3  Attenuation  coefficient  of  the  Si  2p  core  emission  as  a 
function  of  Au  coverage  on  the  Si-Cr  phase.  The  data  show  the  ef¬ 
fect  of  the  presence  of  the  interlayer  on  Si-Au  interdilTusion.  The 
dot-dashed  line  gives  the  corresponding  result  for  the  Si(ll])-Au 
interface  (Ref.  14)  For  Cr  coverage  in  the  monolayer  range- -Hewer 
section),  the  interlayer  slightly  reduces  Si  outdifTusion.  For  Cr  cov¬ 
erages  between  2  and  9  A.  the  interlayer  promotes  Si  outdiffusion. 
while  for  all  Cr  coverages  above  10  monolayers  the  interlayer  sharp¬ 
ly  reduces  the  Si-Au  reaction. 

the  net  effect  of  the  presence  of  the  interlayer  on  the  Si-Au 
reaction  is  done  by  a  normalization  which  takes  into  account 
the  initial  emission  from  Si  surface  atoms  since  this  varies 
with  interlayer  thickness.  Three  different  ranges  of  inter¬ 
layer  thickness  are  clearly  evident.  First,  for  Cr  coverage  of 
0-1  A  or  monolayer  (lower  section  of  Fig.  3).  the  interlayer 
weakly  affects  Si-Au  intermixing  and  reduces  Si  outdif¬ 
fusion.  Secor  J.  for  Cr  coverages  between  1.4  and  9  A.  the 
interlayer  promotes  Si  outdiffusion.  Finally,  for  Cr  coverages 
above  10  monolayers  the  interlayer  sharply  reduces  S i  outdif¬ 
fusion.  The  maximum  promotion  effect  is  obtained  at 
about  2  A  of  interlayer  thickness,  and  the  trend  reverses  it¬ 
self  in  the  monolayer  Cr  coverage  range.  Clearly,  these 
results  cannot  be  interpreted  with  only  a  simple  diffusion 
barrier  model.  Hence,  there  must  be  a  correlation  to 
changes  m  the  local  chemistry  and  energy  balance. 

In  studies  of  SidlP-Cr  interface  formation  we  suggest¬ 
ed1  '  that  there  are  three  formation  stages  for  the  Si-Cr  in¬ 
terface.  In  the  first  («r,  <1-1.5  monolayers),  small  core 
chemical  shifts  and  the  relatively  rapid  attenuation  of  t  ie 
substrate  emission  suggest  little  or  no  interdilTusion  and 
weak  adatom-substrate  interaction  ("weak  chemisorption " ). 
For  !5<4o$.9-l<3  monolayers,  the  core  levels  shift. 


there  is  slower  jticnuation  of  substrate  emission,  and  the 
valence  states  evolve,  tnus  providing  evidence  of  Si-Cr  in- 
terdit'fusion  "reactive  mterdiffusion"1  For  wc.  >  10 
monolayers,  tne  Si-Cr  reaction  is  completed  and  further 
metal  deposition  produces  an  unreacted  Cr  film  on  top  of 
the  reacted  Si-Cr  phase  < "fully  reacted") 

The  clear  one-to-one  correspondence  between  the  dif¬ 
ferent  stages  of  Si-Cr  interface  reaction  and  ine  different  re¬ 
gimes  of  catalytic  effect  of  the  Cr  atoms  in  the  interlayer  is 
compelling  and  allows  us  several  conclusions.  The  promo¬ 
tion  catalytic  effect  of  the  interlayer  on  the  Si-Au  interdif- 
lusion  is  related  to  an  increase  in  reactivity  of  the  Si  surface 
atoms  that  occur  during  the  rejctive  interdiffusion  stage  of 
the  Si-Cr  system.  The  maximum  promotion  effect  occurs 
just  above  the  onset  of  the  reactive  mterdiffusion;  the  inter¬ 
layer  reduces  Si  outdiffusion  below  the  onset  of  the  Si-Cr 
reaction.  The  sharp  transition  suggests  that  during  the  weak 
chemisorption  stage  the  presence  of  Cr  atoms  leaves  the  en¬ 
ergy  content  of  the  Si-Si  bonds  relatively  unchanged  while 
the  metal  atoms  "cover  up"  the  ordered  surface  and  act  as 
a  thin  barrier  against  the  Si-Au  interaction.  When  Cr  and  Si 
begin  io  react,  the  average  Si-Si  binding  energy  is  reduced 
and  the  broken  surface  bonds'*  represent  sites  where  the 
chemically  driven  Si-  Au  intermixing  may  start.  The  chemi¬ 
cally  activated  character  of  this  mterdiffusion  process  is  em¬ 
phasized  by  the  substantial  reduction  of  intermixing  that  oc¬ 
curs  for  interlayer  thickness  above  9  A.  showing  that  even 
very  thin  layers  of  unreacted  Cr  on  top  of  the  Si-Cr  phase 
act  as  an  effective  diffusion  barrier  for  the  intermixing. 

The  gradual  reduction  in  Si  outdiffusion  that  takes  place 
for  increasing  Cr  coverage  between  2  and  9  A  is  related  ei¬ 
ther  to  a  progressive  change  of  the  reactivity  of  the  Si  atoms 
in  tire  Si-Cr  phase  due  to  change  in  bonding  character  or  to 
the  growing  importance,  of  a  rate-limiting  mass  transport 
step  in  the  intermixing  process.  Discrimination  betweerf  the 
two  mechanisms  will  require  temperature-dependent  studies 
(to  vary  the  importance  of  mass  transport  in  limiting  the 
rate  of  reaction)  or  procedure  whereby  the  actual  concentra¬ 
tion  profile  of  Si  atoms  in  the  Si-Cr  interlayer  can  be  deter¬ 
mined.  The  latter  would  allow  a  quantitative  correlation  of 
data  such  as  ours  with  the  effective  surface  concentration  of 
Si  atoms  and  therefore  obtain  the  specific  activity  of  the  Cr 
catalyst  as  a  function  of  interlayer  thickness. 

~  In  summary,  we  have  shown  that  interface  reaction  at  the 
Si-Au  junction  can  be  modulated  by  Cr  atoms  at  the  inter¬ 
face.  that  these  Cr  atoms  can  promote  or  reduce  Si  outdif¬ 
fusion  from  the  bulk,  and  that  the  nonmonotonic  catalytic 
trend  is  related  to  the  existence  of  several  nonequivalent 
microscopic  structures  for  the  Si-Cr  binary  system  at  the  in¬ 
terface.  These  results  point  to  the  importance  of  investiga¬ 
tions  of  the  influence  of  the  interface  catalytic  effect  on  the 
formation  of  the  Schottky  barrier  because  such  studies  could 
clarify  the  relative  importance  of  metal-induced  defects  and 
the  formation  chemistry  of  the  extended  interface  region  in 
determining  the  Schottky  barrier  heights. 
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We  present  synchrotron  radiation  photoemission  studies  of  bulk  CrS^  and  silidde  phases  grown 
on  Si  by  thermal  processing  of  the  Si-Cr  interface.  Experiment  shows  that  Si-Cr  interface  forma¬ 
tion  at  room  temperature  results  in  reacted  phases  that  differ  from  both  bulk  CrSij  and 
in  situ—  grown  Si-rich  CrSij.  Extended— Huckel- theory  linear  combination  of  atomic  orbitals  calcu¬ 


lations  of  the  density  of  states  of  CrjSi,  CrSi,  and 
Si  p  and  Cr  d  states  with  minimal  charge  transfer. 


L  INTRODUCTION 

The  formation  of  refractory  metal  silicides  by  heat 
treatment  of  silicon-metal  interfaces  offers  promise  of 
widespread  application  in  integrated  circuit  technology,1 
but  the  mechanisms  that  drive  interface  formation  and  sil- 
icide  nucleation  are  not  well  understood.2,3  In  part,  this  is 
because  of  limited  microscopic  information  about  the  elec¬ 
tronic  and  structural  properties  of  interfaces  as  they 
evolve  during  thermal  processing.4-7  Further,  there  is 
inadequate  understanding  of  silicides  themselves8  and 
the  parameters  which  influence  their  bulk  properties,  in¬ 
cluding  stoichiometry,  disorder,  segregation,  and  chemical 
bonding.  In  this  paper,  we  focus  on  the  Cr-Si  interface 
and  try  to  resolve  some  of  these  uncertainties  by  examin¬ 
ing  bulk  CrS^,  by  following  the  evolution  of  interface  re¬ 
action  products  as  a  function  of  temperature^  and  by  com¬ 
paring  with  density  of  states  calculations. 

In  an  earlier  investigation  of  the  Si(lll)2X  1-Cr  inter¬ 
face  we  showed  that  room-temperature  reaction  produces 
intermixed  Si-Cr  phases  for  metal  coverages  of  less  than 
10  monolayers.11*12  Electronic  structure  calculations  for 
Cr  silicides  in  simplified  lattice  structure  allowed  us  to 
identify  trends  related  to  bond  formation  and 
stoichiometry  variations.12  The  present  study  extends  that 
work.  Important  new  insight  in  the  silicide  electronic 
structure  is  obtained  through  semiempirical  calculations 
for  CrjSi,  CrSi,  and  CrSi2  in  their  actual  crystal  structure, 
including  the  amount  of  Si-Cr  p-d  hybridization,  the  ener¬ 
gy  position  and  width  of  the  nonbonding  3d  features,  and 
the  role  of  ionicity  in  the  bond.  Experimental  results  for 
bulk  CrSi2  are  then  understood  in  terms  of  substantial  p-d 
hybridization  and  the  reduction  of  the  3d -3d  interaction 
relative  to  pure  Cr.  Comparison  with  the  interface  data 
indicates  both  similarities  and  important  differences. 


CrSij  show  that  Si— Cr  bond  formation  involves 


Temperature-dependent  studies  show  the  evolution  of 
CrSij  and,  ultimately.  Si  segregation  at  the  surface. 

II.  EXPERIMENTAL 

The  bulk  CrSi*  samples  were  prepared  from  ultrapure 
Iochrome  Cr  (99.99%  purity  with  less  than  100  ppm  total 
impurities)  and  Dow  Corning  Si  (resistivity  of  1000 
Q  cm).  The  starting  materials  were  first  arc  melted  then 
directionally  melted  by  if  heating  using  a  horizontal,  cold 
crucible  technique  to  promote  homogeneity  and  grain 
growth.  Analysis  of  the  resulting  CrS»2  ingot  showed 
48.02  wt.  %  chromium  with  lattice  constants  of 
a  =4.4277+0.0001  A  and  c =6.3691  ±0.0004  A  (hexago¬ 
nal  C 40  lattice)  compared  to  the  target  value  of  48.0T" 
wt.  %  and  lattice  constants  of  4.436  and  6.369  A  for  Cr- 
rich-CrSi2  and  4.428  and  6.363  A  for  Si-rich  CrSi2.  (For 
further  crystallographic  information  see  Table  I.) 
Analysis  of  the  microstructure  indicated  small  amounts  of 
Si,  consistent  with  chemical  analysis  and  measured  lattice 
constants.  Gean  CrSi2  surfaces  suitable  for  photoemis¬ 
sion  studies  were  prepared  by  cleaving  posts  of  dimension 
—4x4X10  mm3  in  situ  immediately  before  measure¬ 
ments  were  undertaken. 

The  Si-Cr  interfaces  were  also  prepared  in  situ  by  Cr 
deposition  on  atomically  clean  Sid  11)  surfaces.  Si  wafers 
were  cleaned  for  this  purpose  by  standard  heat  treat¬ 
ments13  and  mild  sputtering.  Interfaces  were  prepared 
with  identical  results  on  both  Si(  1 1 1  )7  x  7  and  Si(  1 1 1 )  1  x  1 
surfaces  as  distinguished  through  photoemission.14  The  in¬ 
terface  behavior  at  room  temperature  was  the  same  as  re¬ 
ported  previously  for  interfaces  on  UHV-cleaved  Si  single 
crystals  in  the  metal-coverage  range  explored  here 
(0>  1— 2  monolayers).11-12  Total  pressure  during  Cr  sub¬ 
limation  from  a  Ta  boat  was  always  better  than  2x  10" 10 
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TABLE  1.  List  of  structural  parameter  for  Cr  silicides. 


Cr2Si 

CrSi 

CrSii 

Crystal  structure 

cubic  A  15 

cubic  520 

hexagonal  C40 

Space  group 

Pm  3m 

P  2,3 

P  6:22 

Lattice  parameters  (A) 

a  =4.5578 

a  =4.607 

a  =4.428;  c=6.363 

Unit-cell  volume  (A5) 

94.68 

97.78 

108.05 

Number  of  molecules/unit  cell 

2 

4 

3 

Number  of  shells  within  7  A 

Si-Si 

3 

9 

20 

Si-Cr 

5 

21 

12 

Cr-Cr 

9 

9 

8 

First-neighbor  distance  (A) 

Si-Si* 

3.947 

2.846 

2.476 

Si-Cr 

2.548 

2.314 

2.476 

Cr-Cr* 

2.279 

2.825 

3.066 

•2.351  A  in  Si. 
"2.494  A  in  Cr. 


Torr  with  recovery  to  the  operating  pressure  of 
~3 X  10-11  Torr  within  ~1  min  after  deposition.  The 
metal  coverage  of  the  Si  substrate,  O,  was  monitored  with 
a  quartz  thickness  monitor  and  is  given  . in  terms  of  the 
Sid  11)  surface  atomic  density  (6=1  monolayer  for 
7.6x  10M  atoms/cm2).  Thermal  processing  of  the  inter¬ 
face  was  performed  by  electron  bombardment  on  the  back 
of  the  Si  wafer.  The  sample  temperature  was  monitored 
with  an  infrared  pyrometer.  Fixed  annealing  temperatures 
were  used  (400  <  T  <  1 100*0,  and  the  annealing  time  was 
increased  in  10- min  cycles  until  the  steady  state-was 
reached,  as  seen  from  the  experimental  spectra.  ■ — 

Photoemission  measurements  were  performed  using 
synchrotron  radiation  from  the  240-MeV  election  storage 
ring  Tantalus  at  the  University  of  Wisconsin- Madison  and 
a  toroidal  grating  monochromator.  Photoelectron  energy 
distribution  curves  (EDCs)  were  recorded  with  a  typical 
overall  resolution  (electrons  plus  photons)  of  0.4  eV  for 
the  valence- band  studies.  The  Si  2 p  and  Cr  3 p  cores  were 
studied  at  hv=  120  eV  (~0.5-eV  resolution)  and  78  eV 
(~  1-eV  resolution),  respectively.  Details  of  the  experi¬ 
mental  system  have  been  given  in  Ref.  15. 

HL  THEORETICAL  TECHNIQUES 

The  energy-band  calculations  for  Cr3Si,  CrSi,  and  CrSi2 
were  based  on  the  linear  combination  of  atomic  orbitals 
(LCAO)  method  in  the  extended— Huckel-theory  approxi¬ 
mation.1*  This  semi  empirical  approach  greatly  simplifies 
the  computation  and  allows  calculations  of  the  electronic 
energy  levels  of  complex  systems  such  as  the  bulk  silicides 
of  interest  here. 17  The  same  calculational  scheme  was  pre¬ 
viously  used  for  the  near-noble-metal  silicides,  yielding  re¬ 
sults  which  were  in  good  agreement  with  experiment. ' 9' 11 
A  detailed  description  of  the  method  can  be  found  in  Ref. 
9. 

In  these  calculations,  the  atomic  valence  orbitals  $a  are 
expressed  as  products  of  the  appropriate  spherical  har¬ 
monics  times  a  radial  wave  function  of  the  form 


R0(r)=c1*a(£l)-|-c,;ra(£2)  , 

where  Ta(£)  are  Slater-type  orbitals  with  orbital  exponent 

The  Slater  orbital  exponents  have  been  adjusted  to 
reproduce  the  electronic  bands  of  pure  Si  (Ref.  19)  and  Cr 
(Ref.  20).  The  resulting  value  of  the  parameters  £,  and  £2 
are  shown  in  Table  II,  together  with  the  constants  Ci  and 
c2.21 

The  Coulomb  integrals  a,~(<^  \H\<j>i)  are  set  equal  to 
the~negative  of  the  valence-orbital  ionization  potential  I, 
for  the  Sf  orbital,22  and  the  resonance  intepals 
Ay  =  <^f  I B  |  d;  )  are  approximated  by9 

Ay  *  —k(I,Ij)i/2Sij  , 

where  Stj  are  overlap  integrals  and  1c =2.5,  as  suggested 
by  Breese  and  Perkins23  for  metallic  systems.  The  J,  can 
be  interpreted  as  charge-dependent  potentials  and  an  itera¬ 
tive  procedure  can  be  built  up  allowing  /,  to  vary  as 

Ii(q)**Ci+Biq  , 

where  q  is  the  excess  charge  per  atom.  This  charge  was 
calculated  through  a  Mulliken  population  analysis,24  and 
the  iterative  procedure  was  stopped  when  the  charge  vari¬ 
ation  between  two  successive  iterations  was  less  than  0.005 
electrons  per  atom.  The  numerical  values  used  for  the  pa¬ 
rameter  C,  and  Bj  (Ref.  25)  for  the  valence  orbitals  are 
also  shown  in  Table  n. 


IV.  EXPERIMENTAL  RESULTS 

Photoelectron  energy  distribution  curves  (EDCs)  for 
the  valence  bands  of  bulk  CrSi2  cleaved  in  situ  are  shown 
in  Fig.  1  for  photon  energies  between  16  and  120  eV.  The 
EDCs  have  been  approximately  normalized  to  the  main 
emission  features  and  are  given  in  arbitrary  units.  The  Si 
Ip  and  Cr  3 p  levels  are  shown  in  the  lower  part  of  the  fig¬ 
ure  for  photon  energies  which  give  comparable  surface 
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FIG.  1.  Valence-band  and  core-level  EDO’s  for  bulk,  polycrystal  CrSij  cleaved  in  situ.  The  measured  binding  energy  of  the  Si  2 p 
(bottom  left)  and  Cr  3 p  (bottom  right)  cores  is  indicated  in  the  figure  and  the  corresponding  experimental  uncertainty  is  0.10  and  0.14 
eV,  respectively. 


sensitivities  for  the  Si  and  Cr  core  levels  (120  and  78  eV). 
The  relative  intensity  of  the  two  main  valence-band  spec¬ 
tral  features  varies  with  photon  energy,  particularly  at  low 
photon  energy  where  final-state  effects  are  important,  but 


both  features  show  very  little  dispersion  for  hv>  25  eV, 
and  we  associate  them  with  structures  in  the  density  of 
states  (DOS).  The  valence-band  emission  of  CrSi2  at 
Av*  50  eV,  where  the  metal  d  photoionization  cross  sec- 
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FIG.  2.  Comparison  of  valence-band  emission  of  CrSi2  at 
hv=  50  eV  with  the  /-projected  partial  DOS  showing  the  Cr- 
denved  id  character.  The  shaded  area  corresponds  to  the  id 
states  that  have  hybridized  with  the  Si  p  states  (in  bonding  com¬ 
binations).  Photoemission  emphasizes  the  nonbonding  d 
features  in  these  3d-metaJ  disilicides.  The  main  spectral  features 
0.6  eV  below  EF  reflects  such  states,  while  important  contribu¬ 
tions  to  the  broad  shoulder  centered  at  1.7  eV  should  come  from 
both  bonding  and  nonbonding  states. 


tion  dominates,  is  compared  with  the  theoretical  Cr- 
derived  3d  character  in  Fig.  2.  The  shaded  region  corre¬ 
sponds  to  d  states  which  are  hybridized  with  Si  p  states. 
(See  the  following  section.) 

The  results  for  the  Sid  1 1)-Cr  interface  are  summarized 
in  Figs.  3  and  4.  In  Fig.  3,  the  valence  band,  Si  2 p  and  Cr 
ip  core  emission  from  the  interface  at  Cr  coverage  0=  10 
monolayers  is  compared  with  the  valence  band  and  core 
emission  of  bulk  CrSb.  Cr  deposition  on  Sid  11)7x7  re¬ 
sults  in  the  formation  of  an  extended  intermixed  Si-Cr 
phase  for  1.5  <0  <  10  which  is  identical  to  what  has  been 
discussed  for  the  Sid  1 1)2  X  1-Cr  interface.  For  0>  10,  a 
film  of  unreacted  Cr  forms  over  the  reacted  interface. 
The  effects  of  annealing  50-A  Cr  films  evaporated  on  the 
Si  wafer  are  shown  in  Fig.  4.  The  topmost  curve 
represents  the  bulk  Cr-like  emission,  and  the  dotted  line 
superimposed  on  this  spectrum  represents  bulk  CrSi2  from 
Fig.  1.  EDC’s  displaced  downward  show  the  effect  of 
10-min  annealing  cycles  at  successively  higher  tempera¬ 
ture  in  the  range  400— 1 100  “C.  Only  at  400  ±50°C  were 
two  heating  cycles  needed  to  reach  steady  state. 


V.  THEORETICAL  RESULTS 

Figure  5  shows  the  total  DOS  for  Cr3Si,  CrSi,  CrSi2  and 
pure  Cr,  as  calculated  with  the  parameters  of  Table  II. 
The  corresponding  DOS  projections  in  Cr-derived  d  and 
s-p  states  are  compared  with  Si-denved  s-p  states  in  Fig. 
6.  These  DOS’s  have  been  calculated  using  a  mesh  of  210 
k  points  in  the  irreducible  Brillouin  zone  for  CrSi2  and  84 
k  'points  for  Cr3Si  and  CrSi.  A  broadening  of  0.1  eV  has 


FIG.  3.  Comparison  of  valence  band  and  core  emission  for  CrSij  to  the  Sid  1 1)-Cr  interface  at  room  temperature  and  Cr  coverage 
0  *  10  ML.  The  photon  energies  have  been  chosen  to  obtain  similar  surface  sensitivity  for  core  and  valence-band  features.  The  inter¬ 
face  reaction  products  (top)  show  a  clearly  different  electronic  structure  with  respect  to  CrSi2,  including  core-binding  energy  differ¬ 
ences  of  —0.76  and  —0.42  eV,  respectively,  for  the  Si  Ip  and  Cr  ip  levels. 
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FIG.  4.  Valence-band  emission  at  /iv=30  eV  of  the  Si-Cr 
junction  as  a  function  of  10-min  annealing  cycles  at  increasing 
temperature.  The  topmost  solid  curve  corresponds  to  a  50-A  Cr 
film  deposited  on  Sill  11)  at  room  temperature.  Superimposed 
on  it  is  an  EDC  for  bulk  CrSi2  (dotted  line).  The  bottom-most 
curve  (dashed  line)  shows  for  comparison  the  theoretical  Cr- 
derived  id  character.  At  — 400*C,  a  disilicide-like  doublet  em¬ 
erges  while  Cr  ip  and  Si  2 p  core  features  show  binding  energies 
consistent,  within  experimental  uncertainty,  with  the  values  ob¬ 
served  for  bulk  CrSij.  Annealing  at  higher  temperatures 
enhances  emission  from  the  Si-metal  p-d  hybrid  states  2—3  eV 
below  EF  (see  Fig.  2)  and  the  Si  Ip  and  Cr  ip  core-binding  ener¬ 
gies  remain  CrSir-like.  The  relative  core  intensity  shows  a  Si- 
rich  stoichiometry  at  the  surface. 


been  used  to  obtain  the  DOS  in  Figs.  5  and  6.  The  excess 
charge  per  atom,9  as  calculated  from  the  iterative  pro¬ 
cedure  outlined  in  Sec.  Ill,  is  the  charge  transfer  per  atom 
due  to  the  Si— Cr  chemical  bond  and  is  given  in  Table  III. 


\s*mj 
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FIG.  S.  Total  densities  of  states  for  Cr  and  Cr  silicides  with 
arbitrary  normalization.  The  corresponding  number  of 
states/(eV  unit  cell)  can  be  obtained  through  the  total  number  N 
of  valence  electrons  per  unit  cell. 


The  total  DOS  for  Cr  is  very  similar  to  that  of  Laurent 
et  al.lb  Strong  id-id  hybridization  splits  the  d  bands  into 
two  groups  of  states,  with  the  bonding  combination  gen¬ 
erally  below  Ef  and  the  antibonding  states  mostly  above 
Ef.  The  Fermi  level  itself  falls  in  a  region  of  low  state 
density.27  Three  main  DOS  features  at  about 
—  1.3,  -2.5,  and  —3.5  eV  correspond  to  states  of  dif¬ 
ferent  symmetry2*  and  spatial  extent.29  Photoemission 
clearly  emphasizes  the  structure  nearest  EF  while  the 
deeper  two  appear  unresolved  at  —  3  eV.30 


TABLE  II.  Input  parameters  used  in  the  calculations  of  the  electronic  structure  of  Cr  silicides.  J 
from  fitting  of  the  energy  bands  of  elemental  Si  (Ref.  19)  and  Cr  (Ref.  20),  C|  and  C:  from  atomic  data 
tables  (Ref.  21),  C,  and  B\  from  Ref.  25. 
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FIG.  6.  Partial  DOS  for  Cr  silkides  showing  the  /-projection 
of  the  Cr-derived  d  (solid  line)  and  s-p  (dashed  line)  character 
versus  the  Si-derived  s-p  character  (dotted  line).  The  chemical 
bonding  depends  on  the  coupling  of  Cr  d  and  Si  p  states.  CrSij 
is  the  only  Cr  silidde  to  have  dominant  nonbonding  DOS 
feature  within  1  eV  of  EF.  The  Si  s  states  are  not  directly  in¬ 
volved  in  the  bonding  and  fall  more  than  8  eV  below  EF . 


VI.  DISCUSSION 

The  chemical  bond  for  the  Cr  and  nearly-noble-metal 
silicides  is  very  similar1- 1011,31  with  coupling  of  Cr  3d 


TABLE  III.  Excess  charge  A?  for  atom  as  calculated  from 
the  iterative  procedure  outlined  in  Sec.  III.  The  procedure  was 
stopped  when  the  charge  variation  between  two  successive  itera¬ 
tions  was  less  than  0.005  electron/atom.  The  charge  transfer  is 
always  less  than  0.5%  of  the  valence  charge/atom  for  Cr  and 
less  than  1  %  for  Si. 


A q  Cr 

A q  Si 

Silidde 

(electrons/atom) 

(dectrons/atom) 

Cr3Si 

-0.013 

+0.038 

CrSi 

-0.029 

+0.029 

CrSij 

+0.015 

-0.007 

and  Si  p  states  and  a  reduction  of  the  Si  s-p  hybridization. 
For  Cr  silicides,  many  of  the  occupied  3d  states  bond  with 
Si,  and  the  resulting  silicide  DOS  is  determined  by  two 
competing  trends.  First,  the  Cr-Si  hybridization  produces 
bonding  features  below  and  antibonding  features  above 
EF,  and  the  3d-band  width  is  broadened.  Second,  the  Cr 
id -Id  interaction  is  reduced,  the  electronic  configuration 
becomes  more  atomic,8,32  and  the  3d-band  width  is  re¬ 
duced.  The  importance  of  the  first  mechanism  in  deter¬ 
mining  the  valence-band  structure  is  greater  than  that  for 
the  nearly-noble-metal  silicides,  because  there  are  few  oc¬ 
cupied  d  states  not  involved  in  the  bonding  with  Si  (non¬ 
bonding  d  states).  For  the  Ni  silicides,  on  the  other  hand, 
the  narrower  nonbonding  3d  states  dominate  the  EDC’s 
and  changes  in  the  Ni-Ni  interaction  determine  most  DOS 
changes  in  the  Ni2Si-NiSi-NiSi2  series.8 

Previous  DOS  calculations  for  Cr  silicides  in  their  teal 
crystal  structures  exist  only  for  the  A 15  structure  Cr3Si 
where  the  total  and  7- projected  density  of  states  as  ob¬ 
tained  by  Arbman  and  Jarlborg33  are  in  good  agreement 
with  those  presented  here.34  In  Cr3Si,  the  Cr-Cr  interac¬ 
tion  dominates  in  determining  the  electronic  structure, 
and  the  total  DOS  in  Fig.  5  is  reminiscent  of  the  3d-3d 
bonding-antibonding  separation  in  bulk  Cr.  Indeed,  the 
Cr-derived  states  between  —4  and  +4  eV  do  not  mix  with 
Si  states.  The  Si— Cr  p—d  bonding  combinations  are  evi¬ 
dent  between  —4  and  —6.5  eV.  The  corresponding  anti¬ 
bonding  states  fall  more  than  4  eV  above  EF.  As  in  other 
metal-rich  silicides,8  the  Si  sp3  hybridization  is  replaced 
by  Si-Cr  p-d  coupling  and  the  s  states  form  an  isolated 
band  around  — 10  eV. 

For  CrSi  (cubic,  £20- type  structure),  the  Si-Cr  p-d 
bonding  states  appear  between  —  1  and  —6  eV,  and  non¬ 
bonding  d  character  is  present  at  the  Fermi  level.  The  p-d 
antibonding  states  are  —4  eV  above  EF,  and  a  gap 
separates  the  s  states  from  the  rest  of  the  DOS,  even  if  the 
value  of  the  energy  separation  is  decreased  with  respect  to 
Cr3Si  because  of  the  larger  Si-Si  interaction.  In  CrSi2 
(hexagonal,  C  40  type)  this  gap  no  longer  exists,  and  states 
more  than  —  7  eV  below  EF  account  for  most  of  the  Si- 
derived  s  character.  The  p-d  bonding  states  appear  be¬ 
tween  —  1  and  —  5  eV,  and  the  empty  antibonding  orbitals 
now  start  only  —  1  eV  above  EF. 

The  results  of  Fig.  5  for  CrSi2  show  that  EF  falls  in  a 
region  of  low,  but  nonzero,  state  density.  Figure  7  shows 
the  DOS  near  EF,  and  the  energy  bands  near  EF  are  repro- 
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FIG.  7.  Detail  of  the  total  DOS  for  CrSij  near  £>.  Top: 
DOS  obtained  with  a  broadening  of  0.1  eV.  Bottom:  High* 
resolution  histogram  plot  showing  the  small  bat  finite  DOS  at 
E,.  — 

dnced  L.  Fig.  8.  The  crossing  of  the  buds  near  M  rales 
out  a  forbidden  gap  and  semiconducting  behavior  for 
CrSi2.  The  carrier  density  in  CiSij  can  be  estimated  by 
counting  the  fraction  of  states  involved  in  the  band  creas¬ 
ing  out  of  the  2038  It  points  in  the  Brillouin  zone.  The  re¬ 
sulting  estimate  of  8. 1  x  1019  carrier/cm3  is  typical  of  a 
semimetaL  (In  our  investigations,  Er  was  determined  ex* 
perimen tally  through  measurements  of  Au  and  Cr  films 
evaporated  in  situ.  The  spectra  of  Fig.  1  are  in  agreement 
with  the  metallic  character  of  CrSi2,  showing  a  low-energy 
cutoff  centered  at  EF  with  a  width  determined  by  the  ex¬ 
perimental  resolution.35  The  low  calculated  carrier  densi¬ 
ty  may  explain  why  CtSi2  has  been  reported  to  be  a  metal 
by  some  authors  and  a  small-gap  semiconductor  by  oth¬ 
ers.3*) 

Our  calculations  show  little  charge  transfer  for  these 
three  riliddes,  always  less  than  1%  of  the  valence 
charge/atom  for  Cr  to  Si  in  Cr3Si  and  CiSi  anc|  complete¬ 
ly  negligible  in  CrSi2  (see  Table  ID.  This  suggests  that 
ioaicity  is  a  minor  contributor  to  the  bond.  Our  calculat¬ 
ed  charge  transfer  for  Cr3Si  appears  in  contrast  to  the 
large  charge  transfer  from  Si  to  metal  suggested  by  Stau- 
demann37  and  Mattbeiss  and  Hamann.31  However,  as  dis¬ 
cussed  by  Pickett  et  a/.,39  numerical  values  depend  on  the 
method  used  to  separate  the  space  around  the  Si  or  metal 


FIG.  8.  Energy-band  diagram  far  CiSij  showing  selected  en¬ 
ergy  levels  near  EF  along  high-symmetry  lines  in  the  C4 0  struc¬ 
ture  Brillouin  zone.  Two  bands  are  shown  to  cross  EF  near  the 
Si  and  L  symmetry  points.  The  resulting  carrier  density  is  es¬ 
timated  at  8. 1 X  101’  carrier/cm3  and  is  typical  of  a  semimetaL 


atoms.  This  problem  does  not  exist  for  our  LCAO 
decomposition  and,  furthermore,  our  calculated  charge 
transfer  is  consistent  in  magnitude  and  sign  with  the  elec¬ 
tronegativity  difference  (Pauling  electronegativity  is  1.6 
for  Cr  and  1.8  for  Si). 

Previous  calculations  far  the  Cr^i,  CrSi,  and  CrSi3  in 
the  CujAu  and  CuAu  structures  showed  qualitatively  the 
same  chemical  bonding  as  the  present  calculations,,  j-e., 
Si-Cr  p-d  coupling,  rtmoval  of  Si  sp 3  hybridization;  and 
small  charge  transfer.12  The  present  results  show  a  gen¬ 
eral  lowering  of  the  W-derived  occupied  DOS  for  CrSi 
and  CrSij.  We  suggest,  that  such  differences  arise  from 
the  use  of  the  real  silidde  geometry  in  the  present  calcula¬ 
tions. 

The  disilidde  of  Cr  is  of  special  interest  since  it  is  the 
only  Cr  silidde  known  to  grow  on  Si  upon  reaction2  at 
45&C.  In  previous  work,  however,  we  demonstrated*  that 
room-temperature  deposition  of  Cr  onto  cleaved  Si  result¬ 
ed  in  intermixing  for  Cr  coverages  up  to  10  monolayers 
(ML).  This  reacted  layer  constitutes  the  “real”  interface 
between  Si  and  Cr  and  is  gradually  covered  with  unreacted 
metal  upon  further  Cr  deposition.  The  implication  is,  of 
course,  that  this  reacted  layer  determines  the  junction  pro¬ 
file  and  transport  properties  of  any  macroscopic  Si-Cr 
junction.  Therefore,  after  examining  in  Fig.  2  the  elec¬ 
tronic  structure  of  CrSi2  on  the  basis  of  our  theoretical 
calculations,  we  compare  in  Fig.  3  the  properties  of  CrSi2 
with  those  of  the  interface  reaction  product  for  6—10 
ML.  The  valence  states  of  CrSi2  show  a  main  peak  at  0.6 
eV  and  a  broad  shoulder  at  about  1.7  eV  (extending  be¬ 
tween  1.3  and  3  eV).  Our  calculations  predict  a  main  non¬ 
bonding  id  feature  0.6S  eV  below  EF,  while  important 
contributions  to  the  broad  shoulder  extending  between  1.3 
and  3  eV  come  from  the  id  bonding  states  at  2.4  eV  in  the 
calculations  and  from  the  nonbonding  id  states  around 
1.8  eV.  The  identification  is  in  agreement  with  the  sys- 
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tematic  trend  in  the  DOS  of  the  3 d  metal  disilicides  (Ref. 
40)  and  confirms  that  ultraviolet  photoemission  em¬ 
phasizes  the  nonbonding  3 d  features  with  respect  to  the  d 
states  directly  involved  in  the  bond  with  Si. 

For  Cr  deposition  onto  Sid  11)  (Fig.  3),  the  Cr-derived 
3 d  states  quickly  dominate  the  valence-band  emission  and 
give  rise  to  a  broad  maximum  that  shifts  with  coverage 
toward  EF  until  the  fully  reacted  layer  is  formed  (6=  10 
ML).  Comparison  of  the  valence  bands  of  the  interface 
reaction  products  and  bulk  CrSi2  shows  that  the  double 
id  structure  (0.6  and  1.7  eV)  characteristic  of  the  disilicide 
is  not  present  for  the  interface  and  only  a  broad  id  peak 
appears  at  0.75  eV  in  the  interface  spectrum  (full  width  at 
half  maximum  of  2.6  eV  compared  to  2. 1  eV  for  CrSi2). 

The  binding  energies  of  the  Si  2 p  and  Cr  ip  cores  are 
important  in  gauging  interface  reaction.  The  Cr  ip  and  Si 
Ip  cores  both  shift  toward  lower-binding  energy  for  cover¬ 
ages  1  <0<  10  and  for  0>  10  stabilize  at  41. 8 ±0. 1  and 
98.80±0. 1  eV  below  EF.4'  Analogous  results  for  CrSi2 
show  a  binding  energy  of  99.54±0. 1  eV  for  the  Si  2 p  and 
42.23±0. 14  eV  for  the  Cr  ip,  identical  within  experimen¬ 
tal  uncertainty  to  the  ip  core  binding  energy  in  bulk  Cr.42 
Therefore,  the  interface- reaction  product  exhibits  Si  Ip 
and  Cr  ip  core  shifts  of  0.74  and  0.42  eV  to  lower  binding 
energy  relative  to  CrSi2,  indicating  that  the  electronic  con¬ 
figuration  differs  substantially  from  CrSi2. 

Examination  of  the  theoretical  Cr-Si  results  of  Fig.  7 
shows  that  only  CrSi2  exhibits  dominant  id  nonbonding 
DOS  features  within  1  eV  of  EF.  Since  the  photoemission 
results  for  all  id  metal  silicides40  are  in  good  agreement 
with  3d-projected  densities  of  states  where  the  nonbonding 
d  character  is  emphasized,  we  conclude  that  none  of  the 
bulk  Cr  silicides  examined  theoretically  exhibits  electronic 
structure  sufficiently  close  to  the  interface-reaction  prod¬ 
uct  to  allow  unambiguous  identification.  Clearly,  what  is 
needed  is  direct  insight  into  the  structural  character  of  the 
interface  so  that  theoretical  modeling  can  be  more  specif¬ 
ic. 

To  examine  whether  the  interface  layer  would  evolve  to 
CrSi2  upon  thermal  treatment,  we  annealed  50-A  Cr  films 
on  Sid  11)  and  monitored  the  valence-band  and  core- level 
features.  In  particular,  annealing  temperatures  in  the 
400— 1100*C  range  were  emphasized  because  Rutherford 
backscattering  and  x-ray-diffraction  investigations2,43 
have  shown  CrSi2  to  be  the  only  growth  phase  to  form  at 
the  interface.  The  topmost  curve  of  Fig.  4  represents  a 
thick  Cr  film  on  Sid  11).  Annealing  at  400±50*C  gives 
strong  Si  Ip  emission  99.45 ±0. 14  eV  below  EF  and  shifts 
the  Cr  ip  cores  from  the  Cr  metal  position  of  42.22  eV  to 
42. 14+14  eV,  converging  on  the  values  of  99.54+0, 10 
and  42.23+0. 14  eV  observed  for  bulk  CrSi2.  Correspond¬ 
ingly,  the  valence-band  emission  becomes  CrSi2-like, 
though  the  spectral  features  are  broader.  Annealing  above 
500*C  results  in  the  emergence  of  enhanced  emission  2-3 
eV  below  EF.  Above  -~700*C,  this  new  feature  dom¬ 
inates,  and  the  id  nonbonding  structure  near  EF  is  only  a 
weak  shoulder. 

These  modifications  of  the  valence  states  result  from 
the  formation  of  a  Si-nch  silicide  at  the  surface. 
Enhanced  emission  2—3  eV  beiow  EF  reveals  an  increased 
number  of  Si-Cr  p-d  hybrid  orbitals  relative  to  the  non¬ 


bonding  d  states.  Further,  the  ratio  of  the  integrated 
emission  from  Si  2 p  and  Cr  ip  cores  is  1.8  times  higher 
for  the  600+  50*C  annealed  interface  silicide  than  for 
cleaved  CrSi2,  indicating  a  marked  enrichment  in  the  Si 
concentration.  Finally,  throughout  the  high-temperature 
annealing  cycles,  during  which  the  valence  states  are 
modified,  the  Si  2 p  and  Cr  ip  cores  maintain  the  bulk 
CrSi2  binding  energy  (experimental  scatter  of  —0.1  eV). 
Low-energy  electron  diffractive  and  Auger  studies  of 
CrSi2  nucleation  on  Sid  11)  by  Oura  er  al.  suggested  that 
annealing  above  600*C— 900*C  produces  several  atomic 
layers  of  silicon  over  a  crystalline  silicide  layer.  Our  re¬ 
sults  confirm  the  increase  in  the  Si  surface  concentration 
but  show  that  the  Si-Cr  bonding  remains  fundamentally 
CrSi2-like,  with  no  evidence  of  formation  of  several  layers 
of  elemental  silicon  at  the  surface.  This  is  also  consistent 
with  analogous  studies  of  Pd2Si.44 

The  results  of  Fig.  6  are  important  in  discussions  of 
room  temperature  Si-Cr  interface  reaction.  We  have 
shown  that  a  CrSi2-like  phase  nucleates  at  relatively  low 
temperature  (400+50*0  and  that  this  phase,  while  lack¬ 
ing  long-range  order,45  exhibits  characteristic  CrSi2-like 
valence  and  core  features.  Therefore,  both  experimental 
and  theoretical  results  indicate  that  the  room-temperature 
interface-reaction  product  is  not  simply  a  disordered  CrSi2 
phase. 

Finally,  we  observe  that  the  Si-rich  CrSi2  phase  shows 
characteristic  enhancement  of  the  valence-band  emission 
2—3  eV  below  EF  while  exhibiting  a  bulk  CrSi2-like  bind¬ 
ing  energy  for  the  Si  2 p  and  Cr  ip  cores.  Therefore,  un¬ 
like  the  case  of  the  Si(  111  )-Pd  interface  where  room- 
temperature  reactions  seem  to  yield  a  locally  Si-rich  Pd2Si 
phase,  we  find  evidence  lhat  room-temperature  Sid  1 1)-Cr 
interface  reaction  does  mot  yield  a  nonstoichiometric-ver- 
sion  of  CrSi2,  the  first  silicide-nucleation  phase. 

The  understanding  of  Si-metal  interface  reaction  re¬ 
quires  the  characterization  of  the  interface-reaction  prod¬ 
ucts  as  a  function  of  the  various  reaction  parameters,  in¬ 
cluding  temperature  and  metal  coverage.  We  have  sought 
to  show  how  a  systematic  experimental  and  theoretical  in¬ 
vestigation  of  the  electronic  properties  of  interface  and 
bulk  Si-metal  compounds  can  increase  our  understanding 
of  both  interface  morphology  and  of  the  parameters,  such 
as  long-range  order  and  stoichiometry  variation,  that  af¬ 
fect  the  silicide-electronic  structure. 
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Synchrotron  radiation  photoemiuion  experiments  with  the  disilicides  of  11,  V,  Nb,  Cr,  Fe,  Co, 
and  Ni  are  combined  with  self-consistent  augmented-spherical-wave  calculations  of  the  density  of 
states  for  metal  siliddes  from  Ca-Si  to  Cu-Si.  These  results  demonstrate  the  importance  of 
silicon p— metal d  bond  formation  extending  to  —6  eV  below  Er  for  all  transition  metals.  Experi¬ 
ment  and  theory  are  combined  to  show  the  movement  of  the  non  bonding  d  states  from  above  Er  for 
Ca-Si  to  well  below  Er  for  Cu-Si.  At  the  same  time,  the  antibonding  Si  p  and  Si  s  states  are  shown 
to  be  relatively  insensitive  to  the  particular  metal  atom  in  the  silidde  series. 


INTRODUCTION  Many  of  the  observations  regarding  the  p-d  bond  are 

therefore  common  to  this  excellent  previous  work.  The 
The  transition-metal  siliddes  have  properties  which  point  of  this  paper  is  to  present  a  systematic  study  of 

make  them  scientifically  interesting  and,  at  the  same  time,  bonding  valid  for  all  of  the  metal  siliddes,  including  the 

extremely  important  for  the  microelectronics  industry.  alkaline-earth  and  near-noble-metal  siliddes. 

Many  of  them  are  chemically  stable  and  resistant  to  corro¬ 
sion  or  degradation,  some  form  at  relatively  low  tempera¬ 
ture  and  have  high  electrical  conductance,1'1  others  show  EXPERIMENTAL  AND  THEORETICAL  TECHNIQUES 
promise  as  interface  diffusion  barriers,1  and  an  increasing 

number  of  siliddes  are  finding  application  in  very-large-  Buttons  of  the  disilicides  of  Ti,  V,  Nb,  Cr,  Fe,  Co,  and 
scale  integrated  (VLSI)  circuit  technology.  Substantial  ef-  Ni  were  initially  prepared  by  arc  melting  high-purity  con- 
forts  are  now  underway  to  understand  the  electronic  prop-  stituents  in  a  nonconsumable  arc  furnace  on  a  water- 
erties  of  siliddes  and  silidde-silicon  junctions.4-7  cooled  copper  hearth.  The  success  and  ease  of  this  process 

We  have  mub-rtulren  a  series  of  experimental- and  varied  from  silidde  to -silidde  but  the  resulting  50-g  sam- 

theoretical  studies  to  delineate  the  role  of  the  Si  s,  tfiTSi  p,  pies  were  generally  solid  and  homogeneous  and  further 

and  the  metal  d  states  in  bonding  and  band-structure  processing  was  limited  to  annealing  to  promote  grain 

properties  of  a  variety  of  transition-metal  siliddes. 1-12  growth.16  In  particular,  this  arc-melting  procedure  was 

This  is  part  of  a  broader  study  which  seeks  insight  not  very  satisfactory  for  TiSi2,  VSij,  NbSi2,  CrSij,  and  CoSi2 

only  about  the  bulk  siliddes  but  also  about  the  electronic  because  they  formed  congruently  from  the  mdt.  The 

and  structural  evolution  of  metal-silicon  interfaces. 12-11  disilidde  NiSi2  forms  peritectically  and  the  reaction  was 

In  this  paper,  we  present  angle-integrated  synchrotron  ra-  completed  by  annealing  at  940*C  for  six  days.  FeSi2  is  ra¬ 
diation  photoemission  results  for  the  bulk  disiliddea  TiSi2l  ported  to  form  congruently  at  1220’C,  but  this  alloy  was 

VSi2,  NbSij,  CrSi*  FeSi2,  CoSij,  and  NiSi2.  Self-  also  heat  treated  for  homogenization  at  940 ‘C  for  13  days, 

consistent  augmented-spherical-wave  (ASW)  calculations  For  CrSi2,  the  initial  arc  melting  was  followed  by  float- 

performed  for  model  metal-silicon  compounds  in  the  zone  melting  to  ensure  uniform  composition, 

high-symmetry  CuAu  and  CujAu  structures  allow  us  to  Samples  of  these  siliddes  were  prepared  for  synchrotron 
model  the  density  of  states  for  disiliddcs  for  comparison  radiation  photoemission  studies  by  fracturing  precut  posts 

with  experiment.  These  results  reveal  the  dependence  of  in  situ  in  the  spectrometer  at  operating  pressures  of 

the  band  structure  on  d-band  occupancy  and  they  demon-  —  5  x  10" 11  Torr.  Radiation  from  the  Tantalus  240-MeV 

strate  the  common  characteristics  of  transition-metal  sili-  electron  storage  ring  was  used  for  studies  in  the  range 

ddes,  including  the  nonbonding  d  states  near  EF,  the  p-d  lO^Av^  140  eV  (Grasshopper  or  toroidal  grating  mono- 

hybridized  states  extending  to  —6  eV  below  EF,  and  the  chromators).  Photoelectron  energy  analysis  was  done  with 

Si  s  states  10—14  eV  below  EF.  For  CoSi2  and  NiSij,  a  commercial  double-pass  electron  energy  analyzer.  The 

which  exhibit  the  CaF2  structure,  results  from  several  dif-  overall  resolution  of  the  experimental  features  varied  from 

ferent  calculations  are  compared  with  experiment.  -0.35  eV,  depending  on  the  photon  energy.  The  results 

It  should  be  pointed  out  that  some  of  the  silicides  dis-  presented  here  are  representative  of  stoichiometric-metal 

cussed  herein  have  been  examined  by  other  authors,  and  disilicides  because  the  samples  were  fractured  in  situ.  Ef- 

Geiatt,  Williams,  and  Moruzzi1  have  reported  detailed  fects  due  to  segregation  at  the  surface  of  silicides  thermal- 

theoretical  studies  of  bonding  between  4 d  transition  met-  ly  grown  on  silicon  or  surface  enrichment  of  one  com- 

als  and  the  non-transition-metal  elements  Li  through  F  portent  due  to  preferential  sputtering  are  thus  avoided.  At 
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the  same  time,  such  effects  are  interesting  because  of  their 
importance  in  interface  phenomena,  and  comparison  with 
our  spectra  is  informative.11,12 

The  experimental  energy  distribution  curves  (EDC’s) 
shown  in  this  paper  have  been  scaled  to  give  an  approxi¬ 
mately  constant  height  for  the  dominant  feature.  At  the 
same  time,  absolute  hv  dependences  of  the  photoioniza¬ 
tion  cross  sections  were  determined  through  normalization 
of  the  emission  intensity  to  the  photon  flux. 

The  phase  diagram  of  silicides  generally  show  a  variety 
of  stable  phases.  In  the  experimental  phase  of  this  study, 
we  emphasize  the  disilicides,  because  the  disilicide  is  the 
most  common  interface  reaction  product  and  because  we 
sought  systematics  in  bonding  within  a  particular 
stoichiometry  range.  The  structural  complexity  of  sili¬ 
cides  has,  no  doubt,  hindered  the  development  of  detailed 
theoretical  pictures  of  their  electronic  structures.17  Few 
calculations  have  been  attempted  which  are,  at  the  same 
time,  descriptive  of  the  correct  stoichiometry  and  crystal 
structure  and  are  also  self-consistent.  The  calculations 
presented  here  are  based  on  the  cubic  CuAu  and  Cu3Au 
structures.  They  show  systematics  in  bonding  as  a  func¬ 
tion  of  stoichiometry  and,  more  important  for  this  paper, 
changes  in  the  Si-derived  and  metal-derived  states  as  the 
metal  itself  is  varied  across  the  3 d  transition  series.  Site 


and  angular  momentum  decomposed  state  densities  de¬ 
rived  from  model  calculations  of  this  kind  have  been 
shown  to  adequately  describe  the  general  trends  of  chemi¬ 
cal  bonding  between  transition  metals  and  non-transition 
metals.  For  example,  the  behavior  of  the  p-d  states  as  a 
function  of  stoichiometry  and  the  estimates  of  charge 
transfer  as  derived  from  such  calculations  and  from  calcu¬ 
lations  using  realistic  silicide  structures  were  found  to  be 
in  remarkable  agreement  for  Ni  and  Pd  silicides.18,19 

The  energy-band  calculations  were  based  on  the 
augmented-spherical-wave  (ASW)  procedure  for  the  solu¬ 
tion  of  the  one-electron  equations.20  Exchange  and  corre¬ 
lation  effects  were  treated  in  the  local-density  approxima¬ 
tion.  Self-consistency  was  pursued  until  the  calculated 
electronic  charges  were  unchanged  to  within  0.001  elec¬ 
trons  within  the  Wigner-Seitz  sphere.  The  calculations 
were  parameter-free,  the  only  a  priori  inputs  being  atomic 
number  and  crystal  structure,  and  they  correspond  to 
theoretical  equilibrium  since  the  atomic  spacing  was  sys¬ 
tematically  varied  until  the  total  energy  was  minimized. 

GENERAL  TRENDS  IN  SILICIDE  BONDING 

In  Fig.  1,  we  show  calculated  densities  of  states  (DOS’s) 
for  the  3 d  transition-metal  monosilicides  based  on  the 


FIG.  1.  Calculated  total  and  partial  state  densities  for  monosilicides  of  Ca  (top)  through  Cu  (bottom)  in  the  CuAu  crystal  struc¬ 
ture.  The  scales  are  given  alonpide  each  curve.  Tic  marks  guide  the  eye  in  observing  systematics  with  increasing  atomic  number,  in¬ 
cluding  the  movement  of  the  d  bands  through  the  Si  p  bands. 
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CuAu  structure.  In  Fig.  2,  we  show  analogous  results  for 
the  trisilicides  (CuAu;  structure).  The  total  DOS  is  shown 
on  the  left  of  each  figure,  the  /-projected  metal  d  character 
is  shown  next,  and  these  are  followed  by  the  projected  sil¬ 
icon  s  and  p  character. 

The  results  of  Figs.  1  and  2  show  how  the  d  bands  dom¬ 
inate  the  DOS  and  move  through  Ep  with  increasing  Z,  as 
indicated  by  tic  marks  at  the  center  of  the  bands.  This 
motion  is  relatively  slow  when  the  d-band  density  of  states 
is  high  at  Ep  (near  the  middle  of  the  transition  series)  but 
is  faster  at  the  extremes  (near  Ca-Si  or  Cu-Si).  The  metal 
d-dehved  states  vary  substantially  within  the  series,  ap¬ 
pearing  first  as  a  broad  band  of  mostly  empty  states  for 
Ca-Si  and  ending  with  the  narrow,  fully  occupied  d  band 
for  Cu-Si.  It  is  important  to  note  that  these  results  always 
show  d  character  well  below  Ef — even  for  Ca  which  has 
s2d°  electronic  configuration  in  the  atom.  This  d  charac¬ 
ter  results  from  the  lowering  in  energy  of  metal-derived  3 d 
states  that  hybridize  with  Si-derived  3 p  states.  The  result¬ 
ing  p-d  bonding  combinations  determine  the  stability  of 
the  silicide,  while  the  corresponding  antibonding  combina¬ 
tions  appear  mostly  above  EF. 

The  Sip— derived  states  are  clustered  into  two  energy 
regions  straddling  the  Fermi  level,  with  the  occupied  or 
bonding  region  being  better  defined  than  the  empty  states. 


As  the  atomic  number  of  the  metal  atom  is  increased  from 
Ca-Si  to  Cu-Si,  the  center  of  the  occupied  cluster  moves  to 
greater  binding  energy  but  the  empty  feature  stays  rela¬ 
tively  fixed  —  3  eV  above  Ep.  This  is  indicated  by  the  uc 
marks  in  the  figures  and  is  best  seen  for  the  trisilicides. 
The  movement  of  the  bonding  states  follows  the  shift  of 
the  metal  d  states  to  greater  binding  energy  with  increas¬ 
ing  Z  and  is  connected  with  the  hybridization  of  silicon  p 
orbitals  with  metal  3 d  states.  The  corresponding  p-d  anti¬ 
bonding  combination  above  Ep  appears  nearly  invariant  in 
energy  in  the  silicide  series. 

The  r-derived  states,  which  extend  to  — 13  eV  below 
Ep,  overlap  relatively  little  with  the  p  states  of  Si  or  the 
metal  d  states.  The  dominant  structure  represents  an  ad¬ 
mixture  of  dominantly  Si  s  and  a  slight  amount  of  metal  s 
character.  For  the  monosilicides,  this  s  band  moves  from 
—  8  eV  for  CaSi  to  — 10  eV  for  CuSi  (see  tic  marks).  For 
the  trisilicides,  it  is  almost  stationary  at  — 10  eV,  except 
for  CaSi]  and  ScSi3.  The  empty  s  states,  on  the  other 
hand,  become  prominent  only  after  Sc-Si  or  Ti-Si,  and  the 
resulting  structure  ultimately  shifts  to  near  Ep. 

The  movement  of  the  d  bands  through  the  Si  p  mani¬ 
fold  in  the  silicide  series  substantially  alters  the  overall 
density  of  states  but  has  relatively  little  impact  on  the  na¬ 
ture  of  the  silicon  p- metal  d  bonds.  The  calculations  indi- 
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FIG.  2.  Calculated  total  and  partial  state  densities  for  the  trisilicidct  of  Ca  (top)  through  Cu  (bottom)  analogous  to  those  of  Fig.  1 
but  with  the  CuAu,  crystal  structure. 
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cate  that  dominant  silicon  p— metai  d  coupling  is  present 
in  all  silicides,  while  the  strength  of  such  coupling  varies 
with  the  relative  position  of  the  p  and  d  states.  In  the 
trisilicide  senes,  for  example,  p-d  bonding  combinations 
appear  in  the  energy  range  0.5— 3.8  eV  for  CaSi3  and  in¬ 
volve  virtually  all  the  calculated  d  charge  per  metal  atom. 
Moving  along  the  id  series,  the  shift  and  sharpening  of 
the  d  bands  yields  p-d  bonding  states  at  1.5—6  eV'in  CrSij 
and  2.5— 6.5  eV  in  CoSi3,  involving  only  about  50%  of  the 
calculated  d  charge  per  Cr  atom  in  CrSi3  and  about  30% 
in  CoSi}.  Correspondingly,  the  metal  d  states  that  are  not 
directly  coupled  with  Si  p  states  form  a  cluster  of  non¬ 
bonding  d  states  nearer  EF. 

We  emphasize  that  the  systematics  developed  through 
Figs.  1  and  2  are  based  on  isostructural  calculations  with 
hybridization  of  certain  states  allowed  by  group  theory 
and  forbidden  to  others.  Calculations  for  different  crystal 
structures  will  modify  the  details  of  the  density  of  states. 
Nevertheless,  we  suggest  that  the  bonding  will  always  in¬ 
volve  p-d  hybridization  and  that  these  hybrid  states  will 
appear  below  the  base  of  the  nonhybridized  states,  in  qual¬ 
itative  agreement  with  our  results.  This  has  been  con¬ 
firmed  by  recent  calculations  for  VSU  (Ref.  21),  CrSU, 
(Refs.  12  and  14),  and  NiSij  (Ref.  19)  in  their  correct  crys¬ 
tal  structures. 
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FIG.  3.  Summary  of  experimental  results  Heft)  and  interpo¬ 
lated  state  densities  for  the  metal  d  (solid  lines)  and  silicon  p 
(dashed  lines)  states  (right)  for  disilicides.  Experiment  verifies 
the  existence  of  metal-silicon  p-d  bonding  states  extending  to  ~  6 
eV  below  Ef,  as  inferred  from  the  calculations,  and  points  to  the 
common  character  of  bonding  in  all  transition-metal  silicides. 
The  dominant  d  emission  near  Ef  corresponds  to  nonbonding 
states.  Note  that  the  results  on  the  right  are  not  scaled  to  true 
state  densities. 


COMPARISON  OF  EXPERIMENT  AND  THEORY 

We  see  from  the  systematics  of  Figs.  1  and  2  that  there 
are  bonding  d  states  below  EF  for  all  these  silicides.  With 
increasing  Z,  the  nonbonding  d  states  sharpen,  move 
through  Ef,  and,  ultimately,  even  the  antibonding  states 
fall  near  EF.  To  facilitate  comparison  of  these  predictions 
with  our  experimental  results,  we  constructed  qualitative 
DOS’s  for  the  disilicides  based  on  results  for  the  monosili- 
cides  and  the  trisilicides.  To  do  this,  we  first  broadened 
the  monosilicides  and  trisilicides  DOS  features  to  retain 
only  major  features.  These  broadened  DOS’s  were  then 
compared  and  a  schematic  DOS  was  drawn  which  had  the 
important  points  common  to  both  but  was,  of  course,  in¬ 
sensitive  to  details.  For  example,  the  d  character  for  CaSi 
was  shown  to  extend  ~4  eV  below  EF  with  a  broad  clus¬ 
ter  of  states  centered  —4.5  eV  above  EF  (width  —6  eV). 


FIG.  4.  Energy  distribution  curves  (EDC’sl  for  TiSij  for 
21  120  eV.  Emission  from  the  Si  2 p  and  Ti  ip  cores  ts 

shown  at  the  top. 
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For  CaSij,  the  d  character  extended  —  5  eV  below  EF  and 
the  empty  d  states  were  centered  —6  eV  above  EF  (width 
3.5—6  eV).  Interpolation  to  CaSi2  would  indicate  d  char¬ 
acter  extending  — 4.S  eV  below  EF  and  empty  states  cen¬ 
tered  —  5  eV  above  EF  (width  ~6  eV).  Analogous  pro¬ 
cedures  were  used  for  the  silicon  p  density  of  states. 
Hence,  on  the  right  of  Fig.  3  we  show  the  schematic 
DOS’s  for  CaSi2  and  for  the  other  systems  under  examina¬ 
tion,  using  dashed  curves  for  the  Si  p  DOS’s  and  solid 
curves  for  the  metal  d  DOS’s.  Note  that  these  DOS’s  are 
sketched  with  arbitrary  units  and  that  the  d  character 
dominates,  as  seen  for  the  quantitative  results  of  Figs.  1 
and  2. 

The  interpolated  p  and  d  DOS's  of  Fig.  3  show  that  the 
d  bands  gather  and  disperse  through  EF  with  increasing  d 
occupancy,  culminating  in  nearly  filled  d  bands  for  NiSi2. 
These  interpolations  predict  clustered  silicon  p  character 
~  5  eV  below  EF.  Hence,  they  preserve  the  most  impor¬ 
tant  aspect  of  the  metal-silicon  bond,  namely,  the  p-d 


Binding  Energy  («V) 


FIG.  5.  EDC’s  for  VSij  for  9^ Av$90  eV  showing  the  hv 
modulation  of  the  initial  sates  of  nonbonding-d  (near  EF),  p-d 
hybrid  (near  4  eV),  p  (near  6  eV),  and  s  (near  10  eV)  character. 


character  well  below  EF.  They  also  establish  the  move¬ 
ment  of  the  nonbonding  cT s  through  EF  with  increasing  Z. 
(These  interpolated  bands  for  NiSi2  predict  a  DOS  d  peak 
— 2.5  eV  below  EF,  compared  to  —2.8  eV  from  the  calcu¬ 
lations  of  Bisi  and  Calandra,19  —  3.5  eV  from  Tersoff  and 
Hamann,22  or  3.3  eV  from  Bylander,  Kleinman,  Rednick, 
and  Grise,23  as  will  be  seen  in  Fig.  10.) 

In  Figs.  4—10,  we  show  EDC’s  for  valence  band  and 
core  emission  for  TiSi2  (Fig.  4),  VSi2  (Fig.  5),8  NbSi2  (Fig. 
6).  CrSi2  (Fig.  7),12  FeSi2  (Fig.  8),  CoSi2  (Fig.  9),  and  NiSi2 
(Fig.  10).9  For  visual  clarity,  the  valence  spectra  have 
been  offset  and  normalized  to  approximately  the  same 
height.  The  EDC’s  for  the  metal  3 p  (4 p  for  Nb)  and  Si  2 p 
core  levels  are  shown  in  the  topmost  or  lowest  sections  of 
the  figures.  Since  most  valence-band  features  in  Figs. 
4—10  do  not  exhibit  major  dependence  on  hv  in  the  range 


Binding  Energy  (eV) 


FIG.  6.  EDC’s  for  NbSij  in  the  photon  energy  range 
10£/iv$90  eV  for  comparison  to  the  3d  transition-metal  disili- 
cides.  Emission  from  the  Nb  4 p  cores  is  shown  at  the  top.  The 
valence-band  spectra  show  r-,  p~,  p-d—  and  d-derived  structures 
analogous  to  those  for  VSi*. 
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10 < hv<  120  eV,  we  associated  them  with  structure  in  the 
density  of  electronic  states. 

The  spectra  in  Fig.  4  for  TiSU  are  dominated  by  intense 
emission  within  3  eV  of  E?  (peak  at  1—1.2  eV)  and  show  a 
shoulder  near  4  eV  for  45  <  h  v  <  65  eV.  At  higher  binding 
energies,  emission  above  the  secondary  background  is  evi¬ 
dent  but  it  is  without  discernible  structure.  (Auger  L  W 
emission  appears  near  9  eV  in  the  40-eV  EDC.)  For  Ti-Si 
most  of  the  d  charge  per  atom  is  directly  involved  in  the 
bonding  with  silicon  and  we  associate  the  main  spectral 
features  to  metal-derived  3 d  states.  The  schematic  DOS 
in  Fig.  3  suggests  that  the  character  of  these  valence  states 
will  become  increasingly  bonding  for  higher  binding  ener¬ 
gies  with  a  parallel  increase  of  the  Si-derived  character  of 
the  valence  states.  The  states  that  account  for  the  DOS 
features  4  eV  below  Ef  exhibit  strong  p-d  hybrid  charac¬ 
ter. 


FlO.  7.  EDCs  for  valence  bands  of  CrSi2  for  1 6  ^  h  v  <;  120 
eV.  Emission  from  the  Si  2 p  and  Cr  3 p  cores  is  shown  in  the 
lowest  part  of  the  figure  (from  Ref.  12). 


The  valence-band  spectra  for  VSU  in  Fig.  5  show  struc¬ 
ture  at  1.6,  5.5,  7.2,  and  10.4  eV.  The  relative  intensity  of 
the  EDC  features  is  modulated  by  the  different  photoioni¬ 
zation  cross  sections  of  initial  states  of  s,  p,  or  d  character. 
On  the  basis  of  the  results  of  Figs.  1  and  2  and  of  Ref.  21, 
we  associate  the  feature  at  1.6  eV  with  metal-derived  3 d 
states,  the  feature  at  5.5  eV  with  hybrid  p-d  orbitJs,  and 
the  structures  at  7.2  and  10.4  eV  with  relatively  pure  Si- 
den  ved  p  and  s  character,  respectively.  These  identifica¬ 
tions  can  be  confirmed  by  a  quantitative  comparison  of 
the  cross  section  for  the  various  initial  states.  Resonant 
photoemission  at  the  Ip  — id  excitation  energy,  in  particu¬ 
lar,  has  already  been  used  for  PdiSi  (Ref.  10)  to  confirm 
the  character  of  the  different  valence  states.  In  the  fol¬ 
lowing  section,  we  will  apply  this  technique  to  VSi2  as  a 
test  case  for  our  systematic  analysis  of  the  disilicide 
valence  states. 


Binding  Enetgy  (eV) 

FIG.  8.  EDCs  for  the  valence  bands  of  FeSi2  for 
21  ^ A 120  eV.  Emission  from  the  Si  2p  and  Fe  3 p  cores  is 
shown  at  the  top. 
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In  Fig.  6,  we  show  EDO’s  for  NbSi2  for  10  <  hv  <90 
eV.  These,  and  results  for  MoSii,*  provide  a  comparison 
for  the  3d  and  4d  isoelectronic  transition-metal  disilicides. 
From  them,  we  see  that  there  are  analogous  spectral 
features  for  NbSi2  and  VSi2  with  slightly  shifted  binding 
energies  (5.2,  8,  and  1 1  eV  for  NbSi2).  The  deepest  can  be 
related  to  the  Si  r-derived  states  and  the  two  at  5.2  and  8 
eV  represent  primarily  p-derived  features  with  substantial 
mixture  of  d  character  in  the  region  of  —4—6  eV.  Fur¬ 
ther,  these  results  show  that  the  d-denved  structure  within 
—4  eV  of  Ef  is  wider  for  NbSii  than  for  VSi2  because  of 
the  inherently  wider  d  bands.  This  d-denved  structure 
closely  resembles  that  of  MoSU  and  the  5 d  transition- 
metal  silicide  TaSi2.* 

Valence-band  and  core  data  for  CrSi2  appear  in  Fig.  7 
(from  Ref.  12).  The  corresponding  calculations  in  Figs. 
1—3  indicate  that  the  d  bands  for  CrSi2  are  approximately 
half  filled  and  that  they  draw  away  from  EF  as  they  be¬ 
come  nearly  filled.  Indeed,  for  NiSi2  the  d  bands  are 
predicted  to  be  —  3—6  eV  below  Ef.  The  effect  of  the  d- 
band  filling  can  be  seen  by  examining  the  CrSi2,  FeSi2, 
CoSi2,  and  NiSi2  silicide  senes  (Figs.  3  and  7 — 10).  The 
EDCs  for  CrSi2  indicate  strong  emission  near  EF  (0.6  eV), 
structure  near  —2  eV,  and  a  shoulder  at  —6  eV.  The 
EDO’s  for  FeSi2  show  less  structure  (peak  at  about  0.75 
eV  and  broad,  structureless  emission  extending  to  about 
- 10  eV).  EDO’s  for  OoSi2,  on  the  other  hand,  show  the 
dominant  d-derived  feature  to  have  drawn  away  from  EF, 
as  predicted  by  the  interpolated  DOS  results  of  Fig.  3 
(peak  at  1.4— 1.5  eV  and  shoulder  at  —3.8  eV  below  EF). 
Results  for  NiSi2,  shown  in  Fig.  10  for  hv= 40  and  50  eV, 
are  very  similar  to  those  for  OoSi2  except  that  the  location 


Binding  Energy  (eV) 

FIG.  9.  EDCs  for  the  valence  bands  of  CoSij  in  the  photon 
energy  range  25  120  eV  compared  to  the  total  density  of 

states  calculated  by  Tenoff  and  Hamann  (Ref.  22). 


of  EF  is  higher  within  the  valence  bands,  the  d  bands  are 
almost  completely  filled,  and  the  separation  of  the  non 
bonding  and  bonding  states  is  reduced. 

These  various  results  can  be  understood  on  the  basis  ol 
the  systematics  of  the  schematic  DOS  of  Fig.  3,  taking 
into  account  the  empirical  observation  that  photoemission 
spectra  for  10  <  /t  v  <  120  eV  emphasize  nonbonding  3d 
states  over  the  hybrid  p-d  orbitals.  This  can  be  seen,  for 
example,  by  comparing  the  NiSi2  calculations  of  Fig.  10 
with  experimental  results  or  comparing  our  synchrotron 
radiation  results9  with  XPS  (x-ray  photoelectron  spectros¬ 
copy)  results  by  Grunthaner,  Madjukar.  and  Mayer  for 
Ni2Si.7  In  our  EDO’s  the  DOS  structures  3—3.5  eV  below 
EF  for  NiSi2  (2  eV  for  CoSi2)  dominate  the  spectra  but  the 
5-eV  feature  (4  eV  for  CoSi2)  is  much  weaker  than  in  the 
XPS  results  and  in  the  calculations.  Bearing  this  in  mind, 
the  experimental  emission  feature  at  0.6  eV  for  CrSi2  can 
be  associated  with  nonbonding  metal  3d  states,  while 
emission  at  2.0  eV  and  below  should  come  from  bonding 
p-d  orbitals.  The  shoulder  at  6  eV  corresponds  to  DOS 
features  of  substantial  Si  —p  character.  Similarly,  for 
FeSi2>  where  a  large  fraction  of  the  metal  3d  charge  is  not 
coupled  with  Si-derived  states,  the  main  emission  features 
within  -2  eV  of  EF  reflect  nonbonding  metal  d  states, 
while  the  broad  shoulder  —  5  eV  below  EF  exhibits  mixed 
p-d  character. 

CoSi2  is  near  the  end  of  the  3d  row  where  the  nonbond- 


FIG.  10.  EDO's  for  the  valence  bands  of  NiSij  for  )i  v=40 
and  50  eV  (from  Ref.  9).  For  comparison  we  show  that  total 
density  of  states  calculated  by  Bisi  and  Calandra  (Ref.  19)  (top), 
Tersoff  and  Hamann  (Ref.  22),  and  Bylander,  Kleinman,  Med- 
nick,  and  Grise  (Ref.  23)  (bottom). 
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mg  d  states  form  a  sharp,  dominant  emission  feature  at 
1.45  eV  and  the  bonding  d  states  appear  as  a  relatively 
weak  structure  at  about  2-eV  higher  binding  energy.  In 
Fig.  9  we  show  the  experimental  results  for  CoSi,  and 
compare  them  to  the  calculated  density  of  states  of  Ters- 
off  and  Hamann.22  Inspection  shows  quite  good  agree¬ 
ment  but  indicates  that  the  calculations  overestimate  the 
binding  energy  of  the  experimental  feature  by  0.3— 0.5  eV. 

NiSi2  is  probably  the  most  studied  silicide  both  from 
the  experimental  and  theoretical  points  of  view.  Angle- 
integrated  and  angle-resolved  photoemission  investigations 
of  its  electronic  structure  have  recently  been  reported  by 
ourselves,9  Chang  and  Erskine,  and  Chabal,  Hamann, 
Rowe,  and  Schliiter.24  In  Fig.  10  we  summarize  the 
angle-integrated  results9  and  compare  them  to  the  densi¬ 
ties  of  states  calculated  by  Tersoff  and  Hamann,22  Bisi 
and  Calandra,19  and  Bylander,  Kleinman,  Mednick,  and 
Grise.23  As  shown,  there  is  good  overall  quantitative 
agreement  with  experiment  and  the  different  sets  of  calcu¬ 
lations  yield  very  similar  densities  of  electronic  states,  al¬ 
though  the  details  of  the  bonding  vary  from  author  to  au¬ 
thor  (charge  transfer  and  ionic  contribution  to  the  bond, 
electronic  configuration  at  equilibrium  for  the  metal  and 
Si  atoms). 

RESONANT  PHOTOEMISSION 

Identification  of  ,he  character  of  the  valence  states  can 
be  supported  by  quantitative  analysis  of  the  cross  sections 
of  the  various  initial  states.  To  do  this,  we  show  in  Fig.  5 


FIG.  11.  Photoionization  cross  sections  for  VSi;  correspond¬ 
ing  to  initial  states  marked  by  arrows  in  Fig.  5.  As  a  result  of 
3p-3d  resonance,  the  spectra  show  a  threefold  enhancement  of 
the  nonbonding  d  states.  For  the  p-d  hybrid  states,  the  magni¬ 
tude  of  the  enhancement  depends  on  the  degree  of  hybridization 
and  decreases  for  increasing  Si  p  contribution. 


a  series  of  EDO’s  for  VSi2  and  the  results  of  analysis  of 
the  peak  heights  as  a  function  of  photon  energy  in  Fig.  11. 
VSi2  was  chosen  because  several  structures  could  be 
resolved  in  the  valence  bands,  including  features  at  1.6  eV 
(a  derived),  5.5  eV  (p-d),  7.2  eV  (p),  and  10.4  eV  (s  de¬ 
rived).  Emission  from  the  Si-derived  s  and  p  states  is  low 
for  hv— 50— 60  eV  but  then  increases  at  higher  energy. 
The  d  bands,  on  the  other  hand,  are  always  prominent,  but 
this  is  partially  an  artifact  of  normalization  of  the  spectra 
of  Fig.  5. 

Partial  photoionization  cross  sections  were  determined 
from  the  spectra  of  Fig.  5  and  others  for  9  < /tv  <90  eV 
measured  at  increments  of  1  or  2  eV.  These  are  shown  in 
Fig.  1 1  for  initial  states  at  0.3,  1.6,  3.2,  5.5,  and  7.4  eV  (see 
arrows  in  Fig.  5).  Gear  enhancement  of  these  states  (and 
the  entire  band)  for  hv>  38  eV  is  associated  with  the  exci¬ 
tation  of  the  V  Ip  core  (Eg— 37 A  eV).  For  E,  =0.3  and 
1.6  eV,  the  emission  decreases  for  20<Av<40  eV  but  is 
then  enhanced  approximately  threefold  between  -  38  and 
—  50  eV.  For  states  in  the  p-d  bonding  region,  the 
enhancement  is  weaker  '  it  the  shape  is  approximately  the 
same  for  hv>  38  eV.  States  with  increased  Si-derived 
character  exhibit  considerably  less  enhancement.  (The 
structure  that  appears  around  hv=40  eV  is  due  to  over¬ 
lapping  Auger  emission.) 

Resonant  valence-band  enhancement  has  been  observed 
in  other  silicide  systems10  and  is  analogous  to  d-f 
enhancement  studied  in  lanthanides  and  actinides.25  This 
resonance  reflects  the  quantum-mechanical  equivalence  of 
different  paths  leading  from  the  ground  state  to  the  final 
state  where  an  emitted  electron  is  in  the  continuum.  The 
two  paths  are  direct  valence-band  emission 
3p6(4sd)"  +  kv-*3p<,(4sd)"-1 +e  and  a  core-lcvel- 
involving  process  3pb(dsd)*  +  hv—3pi(dsd )*  ^ 1 
— •  3pt,{Asd)*~'  +e.  The  strength  of  the  enhancement  is  a 
measure  of  the  overlap  of  the  p- core  hole  and  the  d- 
valence  state.  Hence,  as  the  d  states  become  hybridized 
the  overlap  is  reduced  and  the  enhancement  diminishes. 
The  importance  of  this  phenomenon  to  us  here  is  that  it 
can  be  used  to  highlight  d  emission  in  the  photon  energy 
range  near  the  p-core-level  excitation  energy.  The  results 
for  VSi2  (as  for  Pd2Si  in  Ref.  10)  clearly  show  the  dif¬ 
ferent  character  of  Si-  and  metal-derived  states  and  they 
support  the  identification  of  the  spectral  feature  as  3d  de¬ 
rived,  with  mainly  nonbonding  character  at  EF  and  in¬ 
creasing  p-d  hybrid  character  for  increasing  binding  ener¬ 
gy  in  the  energy  range  2—6  eV. 

Cross-section  analyses  for  MoSi2  similar  to  those  for 
VSi2  show  similar,  but  smaller,  enhancement  for  4 p-Ad  in¬ 
teraction  as  for  3p-3d.  Analogous  studies  for  ScH2,  YH2, 
LaH2,  and  ThH2  (Ref.  26)  indicate  reduced  enhancement 
for  increasing  principal  quantum  number  n  in  np-nd  reso¬ 
nance.  This  is  due,  in  part,  to  the  increasing  delocaliza¬ 
tion  of  the  d  wave  function. 

CONCLUSIONS 

In  this  paper  we  have  given  experimental  and  theoreti¬ 
cal  insight  into  bonding  systematics  for  the  metal  disiii- 
cides  over  the  full  range  of  the  3d  transition  metals.  The 
systematics  indicate  the  following: 
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(1)  The  states  which  we  see  as  dominant  in  synchrotron 
radiation  photoemission  are  generally  nonbonding  id 
states  since  they  overlap  very  little  with  Si-derived  states. 
The  results  for  the  disilicides  of  Ti  through  Co  all  show 
important  d  character  within  3  eV  of  EF  (lower  for  Ni). 
Hence,  the  electrical  properties  of  silicides  will  be  deter¬ 
mined  by  d  states. 

(2)  Silicon/?— metal  d  mixing  is  important  for  all  sili¬ 
cides  and  lowers  metal  d  states  to  well  below  EF.  Our  ex¬ 
perimental  observation  of  bonding  mixtures  of  p-d  states 
supports  the  model  that  p-d  bonding  states  will  generally 
account  for  the  silicide  stability.  Detailed  calculations  of 
the  future  must  consider  the  specifics  of  the  p-d 
hybridization — this  must  vary  with  stoichiometry,  crystal 
structure,  and  possibly  nonequivalent  sites  (see,  for  exam¬ 
ple,  the  calculation  of  Switendick27  for  Fe3Si). 

(3)  The  detection  of  the  Si  s-derived  states  has,  so  far, 
been  limited  to  the  disilicides  of  V,  Nb,  Mo,  Ta,  and  Ce. 
(Parks,  Reihl,  Martensson,  and  Steglich28  have  observed 
the  onset  of  the  Si  s  band  near  12  eV  for  CePd2Si2  and 
CeAu2Si2.)  For  the  other  disilicides,  we  observe  only  an 


unstructured  tailing  off  of  valence-band  emission  near 
10—12  eV.  The  invariance  of  this  feature  indicates  that 
the  Si  s  character  is  relatively  unimportant  in  bonding. 

We  anticipate  that  these  results  can  be  generally  inter¬ 
polated  to  the  monosilicides  [e.g.,  MnSi  and  NiSi  (Ref.  9)] 
and  the  metal-rich  silicides  [Fe3Si,  Ni2Si  (Ref.  9),  Pd2Si 
(Ref.  10)].  They  will  most  likely  find  substantial  validity 
for  the  germanides  as  well. 
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Cluster-Induced  Reactions  at  a  Metal -Semiconductor  Interface:  Ce/Si ( 111 ) 
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Abstract 

Synchrotron  radiation  photoemission,  angle-resolved  Auger,  and  LEED  studies 
show  Ce  cluster  formation  on  Si (111),  These  nonmetal  lie  clusters  grow  for 
Coverages  0.1  to  0.6  ml,  interact  weakly  with  the  substrate,  and  induce 
200  meV  band  bending  changes.  At  ~0.6  ml,  they  stimulate  surface 
disruption,  producing  a  metallic  interfacial  silicide.  The  association  of 
d-  or  f-band  metal  clusters  with  surface  reaction  substantially  extends  the 
cluster-induced-reaction  model  proposed  for  Al/GaAs. 


PACS:  73.40.+y 

73.40.  -c 

73. 40.  Ns 
79.60.Gs 


-2- 


Microscopic  examinations  of  evolving  metal /semi conductor  interfaces 
have  revealed  several  systems  that  react  only  when  the  metal  coverage,  0, 
exceeds  a  critical  coverage,  ec.1-10  Modeling  of  these  interfaces  requires 
an  understanding  of  the  triggering  phenomenon.  Despite  considerable 
efforts,  however,  the  mechanism  remains  elusive. 

Zunger1*  has  suggested  that  cluster  formation  can  induce  defects  at  the 
Al/GaAs(110)  interface.  His  total  energy  calculations  stressed 
adatom-adatom  interaction  through  clustering,  a  point  of  view  which  differs 
from  those  which  stress  metal /substrate  bonding.1^*15  In  light  of  the 
significance  of  this  cluster  mechanism,  we  sought  a  metal /semi conductor 
system  in  which  the  onset  of  chemical  reaction  could  be  linked  to  the 
ripening  of  clusters. 

In  this  paper,  we  describe  the  room  temperature  evolution  of  Ce 
clusters  on  Si ( 111 ) -2x1  and  7x7  surfaces.  High  resolution  synchrotron 
radiation  photoemission,  angle-resolved  Auger  spectroscopy  (ARAS),16  and 
LEED  measurements  show  changes  in  both  electronic  structure  and  morphology. 
Indeed,  this  combination  of  techniques  should  make  it  possible  to  assess 
cluster  formation  for  other  interfaces  and  should  reveal  the  relative 
importance  of  adatom/substrate  and  adatom/ adatom  bonding. 

The  photoemission  experiments  were  performed  at  the  Wisconsin 
Synchrotron  Radiation  Center  using  "Grasshopper"  and  toroidal  grating 
monochromators  for  12<hv<135  eV.  Photoelectrons  were  analyzed  by  a  double 
pass  CMA  with  overall  energy  resolution  of  between  200  meV  (hv  *  40  eV)  and 
500  meV  (135  eV).  -The  ARAS  studies  used  a  single  pass  CMA  with  limited 
electron  acceptance  and  a  precision  manipulator  to  vary  the  take-off  angle 
from  normal  emission  to  grazing  incidence.  This  detector  was  also  used  to 
obtain  I-V  LEED  profiles  along  high  symmetry  lines. 


Clean  Si ( 111 ) -2x1  and  7x7  surfaces  were  produced  by  cleavage  or 
standard  A r+  sputtering  and  annealing.  High  purity  cerium  was  evaporated  at 
pressures  better  than  lxlO"18  Torr;  the  thickness  was  monitored  by  an 
oscillating  quartz  crystal  and  the  evaporation  rate  was  ~lA/min.  Operating 
pressures  in  both  spectrometers  were  3-5xl0-11  Torr. 

To  demonstrate  the  existence  of  a  critical  coverage  corresponding  to 
the  onset  of  reaction,  we  show  in  Fig.  1  Si  2p  core  level  spectra  taken  at 
hv*135  eV  (escape  depth  ~4A).17  For  cle^n  SI  ( 111 )-2xl ,  the  2p  doublet  is 
broadened  by  the  surface-shifted  component.17  Between  zero  and  -0.6  ml, 
there  is  a  rigid  shift  of  200  meV  due  to  variations  in  band  bending  (as 
confirmed  by  bulk  sensitive  spectra  at  hv*106  eV),  but  the  unvarying 
lineshape  indicates  minimal  Ce-substrate  interaction.  Analysis  of  the  Si-2p 
lineshape  for  e>0.6  ml  reveals  a  reacted  Si  species  shifted  0.7  eV  to  lower 
binding  energy.  This  component  grows  and  is  most  pronounced  near  3  ml  and 
the  substrate  component  attenuates  rapidly  (discussed  elsewhere  using 
integrated  intensities  of  the  different  components18).  From  coverage 
studies  with  0.06  ml  increments  between  0  and  1  ml ,  we  conclude  that 
ec*0.6  ml . 

LEED  studies  showed  conversion  from  2x1  or  7x7  patterns  to  lxl  by 
-0.4  ml  but  showed  no  metal -induced  superstructure.  The  Si -lxl  pattern 
persisted  through  ec,  although  with  increasing  background.  For  0.8  ml, 
where  photoemission  showed  intermixing,  LEED  showed  sufficient 
uncovered/unreacted  Si  for  coherent  diffraction  (lxl  pattern).  This 
indicates  that  the  reaction  could  not  be  proceeding  uniformly  on  the 
surface,  and  points  to  cluster-induced  reaction.  Even  at  1.2  ml  a  lxl 
pattern,  which  was  scarcely  visible  to  the  eye,  could  be  detected  by  taking 
I-V  measurements  in  high-symmetry  azimuthal  planes  perpendicular  to  the 


-4- 


surface.  By  2  ml ,  the  (10),  (11),  and  (20)  LEED  beams  were  completely  gone, 
indicating  the  total  destruction  of  long  range  order. 

Valence  band  studies  at  hv*35  and  60  eV  were  undertaken  to  examine  the 
electronic  states  at  low  coverage.  Changes  induced  by  d-derived  states  were 
readily  seen  at  35  eV;  comparison  of  results  for  hv*35  and  60  eV  showed  the 
emergence  of  Ce  4f  emission. *®*2°  Shown  in  Fig.  2,  the  Ce-derived 
emission  at  -3.3  and  -1.2  eV  for  0.1<e<0.6  ml  rapidly  overwhelms  Si-emission 
and,  with  increasing  coverage,  the  feature  at  -1.2  eV  grows  in  relative 
intensity.  Simultaneously,  the  FWHM  measured  at  hv*35  eV  diminishes  from 
-4.4  eV  for  0*0.2  ml  to  3.9  eV  for  0*0.4  ml.  Consistent  with  the  core  level 
results,  we  find  that  at  0*0.6  ml  the  valence  band  center  shifts  abruptly 
toward  Ep  by  -0.7  eV  and  broadens  sufficiently  to  lose  its  prominent  doublet 
character.  The  appearance  of  the  Fermi -level  cutoff  indicates  metallic 
character  for  the  reacted  phase. 

The  valence  band  results  indicate  cluster  formation  and  growth  for 
coverages  0-0.6  ml.  Adatom-induced  states  are  clearly  visible  but  are 
shifted  away  from  Ep,  analogously  to  what  has  been  observed  in  all 

pi 

photoemission  studies  of  weakly  interacting  clusters  on  inert  substrates. 

The  observed  cluster  bandwidth,  which  is  substantially  greater  than  that  of 
bulk  Ce,*®  suggests  a  large  number  of  inequivalent  atoms,  as  observed  for 
metal  clusters.  The  absence  of  a  Fermi  level  for  0*0.6  ml  indicates  that 
the  molecular  cluster  has  not  reached  the  metallic  droplet  limit. 

Studies  of  Ce  4f  emission  probe  the  environment  around  the  Ce  atom  as 
effectively  as  Si  2p  emission  reveals  the  Si  environment.  Modification  of 
the  Ce  environment  above  ec  can  be  seen  through  difference  spectra  obtained 
by  subtracting  valence  band  EDCs  at  35  and  60  eV  (Fig.  2).  This  procedure 
emphasizes  the  4f  character,  as  shown  for  a-Ce,  y-Ce,*®  and  the  Ce 
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pnlctides.  Even  for  Ce  coverages  of  0.1  ml,  a  clearly-defined  spectral 
feature  can  be  associated  with  the  Ce  4f  at  -3.3  eV.  The  changing  4f 
fingerprint  can  therefore  be  used  to  assess  chemical  reaction  (4f  final 
states  screening  will  be  discussed  elsewhere). 

Figure  2  shows  that  the  4f  feature  sharpens  but  does  not  shift  with 
coverage  below  ec.  At  e*0.8  ml  its  appearance  is  modified  by  a  shoulder  at 
0.7  eV  lower  binding  energy,  and  by  1.2  ml  the  4f  feature  has  fully  shifted 
to  the  lower  binding  energy  position.  Comparison  with  results  obtained  from 
a  bulk  sample  of  CeSi2  shows  a  silicide  formation. 

At  the  top  of  Fig.  3  we  show  Ce(N^VV)  and  Si^^VV)  integrated  Auger 
peak  intensities  measured  at  normal  emission  as  a  function  of  coverage. 

With  increasing  coverage,  the  Ce  signal  steadily  increases  while  the  Si 
signal  is  attenuated  significantly  from  0  to  0.4  ml.  Between  0.4  and 
0.8  ml,  corresponding  to  the  transition  through  ec,  the  SI  attenuation  is 
reduced.  This  behavior,  which  deviates  from  that  expected  for  either 
ordered  overlayers  or  an  intermixed  phase,  can  be  explained  by  initial  Ce 
coverage  of  a  fraction  of  the  substrate,  followed  by  an  intermediate  region 
0.4-0. 8  ml,  where  two  processes  compete.  Initially,  cluster  formation 
attenuates  the  substrate,  but  this  is  offset  by  Si  intermixing  once  reaction 
starts.  The  reacted  regions  then  spread  laterally,  consuming  more  of  the 
surface  and  attenuating  the  Si  emission. 

Polar-angle  profiles  of  the  Ce  and  Si  Auger  Intensities  measured  in  the 
(110)  azimuthal  plane  are  shown  in  the  lower  two  sections  of  Fig.  3. 
Structure  in  these  polar  profiles  would  reflect  electrons  scattering  from 
near  neighbor  atoms  and  subsequent  interference.  For  Ce,  the  results  are 
featureless,  consistent  with  disordered  clusters.  In  contrast,  the  polar 
profile  for  clean  Si  shows  diffraction  modulation  characteristic  of  the  7x7 
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surface,  and  this  modulation  remains  in  detail  at  0.8  ml.  Since  atomic 
scattering  factors  for  Ce  are  roughly  three  times  thos*  for  Si  for  100  eV 
electrons,  a  uniform  Ce  layer  should  obscure  fine  structure  related  to 
substrate  emission.  The  persistence  of  substrate  structure  indicates  that  a 
substantial  fraction  of  the  Si  electrons  are  unscattered,  i.e.  much  of  the 
surface  is  free  of  Ce.  (Results  for  intermixed  Au/Si  show  Si  Auger 
diffraction  fine  structure  to  be  damped  substantially  more.)  On  the  other 
hand,  the  Si  attenuation  reflects  inelastic  scattering  of  electrons  which 
propagated  through  the  Ce  clusters.  The  damping  of  the  fine  structure  at 
1.2  ml  indicates  the  increasing  size  of  the  reacted  regions,  consistent  with 
cluster  induced  reaction.  By  coverages  of  4  ml ,  the  Si  polar  profile  is 
completely  featureless. 

We  have  shown  here  that  clusters  form  by  ~0.1  ml  coverage,  but  produce 
Ce/Si  intermixing  only  for  3>0.6  ml.  Hence,  the  mechanism  proposed  by 
Zunger  for  Al/GaAs  is  not  limited  to  simple  metals  on  compound 
semiconductors,  but  may  be  important  for  a  broad  class  of  materials.  We 
hope  that  this  work  will  stimulate  theoretical  efforts  to  describe  the  early 
stages  of  interface  formation,  including  clusters  where  appropriate,  and 
establish  the  theoretical  basis  for  the  metal -overlayer-triggered  reaction. 

This  work  was  supported  by  ARO  DAAG29-83-K-0061  and  a  Northwest  Area 
Foundation  Grant  of  the  Research  Corporation  (S.A.  Chambers).  Discussions 
with  G.  Margaritondo,  F.  Grunthaner,  and  A.  Franciosi  are  gratefully 
acknowledged. 
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FIGURE  CAPTIONS 


Fig.  I  Si  2p  core  total  emission  and  deconvolution.  The  broken  vertical 
line  denotes  200  meV  change  in  band  bending.  Representative 
statistics  are  given  by  dots  for  the  clean  surface  at  bottom. 

Fig.  2  Energy  distribution  curves  for  hv*35  and  60  eV  (left)  and  the 

difference  curves  (right)  to  highlight  the  Ce  4f  character.  Below 
0.6  ml,  the  4f  indicates  cluster  formation.  Above  0.6  ml,  the  4f 
features  shifts  0.7  eV  and  reflects  local  coordination  with  Si,  as 
shown  by  comparison  with  CeSi2  (bottom). 

Fig.  3  Attenuation  curves  for  Ce(N4VV)  and  Si^^VV)  Auger  intensity 

from  normal  emission  spectra  (top.  Polar  profiles  are  shown  for  the 
Ce(N4VV)  and  Si (L2  3VV)  Auger  intensities  for  different  coverages 
(middle  and  bottom).  90°  is  normal  emission.  All  spectra  taken 
in  the  (110)  azimuthal  plane  perpendicular  to  the  surface. 
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QUANTITATIVE  INTERDIFFUSION  STUDIES  OF  NOBLE  METAL/Si (111 )-7x7  INTERFACES 
BY  ANGLE -RESOLVED  AUGER  ELECTRON  EMISSION 

I 
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ABSTRACT 

Polar-angle  resolved  Auger  electron  emission  and  related  modeling  have 
been  used  to  determine  quantitative  details  of  composition  and  diffusion 
phenomena  in  the  Cu/Si ( 111 )  and  Au/Si (111)  interface  systems.  The  extent  of 
Si  outdiffusion  increases  with  coverage  in  a  gradual  fashion  for  Au,  with 
characteristic  diffusion  lengths  (CDL)  ranging  from  18%  of  the  overlayer 
thickness  at  a  2A  coverage  to  39%  at  20A.  In  contrast.  Si  CDL  values  in  Cu 
increase  abruptly  from  20%  of  the  overlayer  thickness  at  2A  to  -35%  at 
coverages  from  5  to  20A.  Au  diffuses  deeply  into  the  Si  substrate,  whereas 
Cu  indiffuses  only  slightly  (CDL  =  1-2A).  These  model  predictions  are 
consistent  with  the  more  substantial  lattice  disruption  brought  about  by  Au, 
leading  to  no  long-range  order  in  the  overlayer,  and  the  persistent 
diffusion  of  Si  to  the  surface  of  150A  thick  Au  overlayers. 

PACS:  64.75.+g 
66.30. Ny 
73.40. -c 
79.60.Gs 
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Evaporation  of  metal  overlayers  on  single-crystal  semiconductor 
surfaces  at  room  temperature  often  results  in  substantial  interdiffusion,  a 
phenomenon  not  frequently  observed  in  metal /metal  or  gas/metal  chemisorption 
systems.  Quantitative  characterization  of  such  interfaces  places  a  new 
challenge  before  experimental ists.  Although  ion  sputter  depth  profiling  may 
be  useful  for  semi -quantitative  characterization  of  thick  overlayers,  its 
utility  in  studying  thin  overlayers  (<20A)  is  limited.  It  seems  that  a  non¬ 
destructive  method  of  depth  profiling  that  is  sensitive  to  ultra-thin  metal 
overlayers  would  be  of  considerable  value.  In  this  article  we  describe 
the  first  application  of  polar-angle  resolved  Auger  electron  emission  to  the 
quantitative  characterization  of  metal -Si  interfaces.  We  have  chosen  the 
noble  metals  Cu  and  Au  because  both  are  qualitatively  known  to  Intermix  with 
Si. i-4 

We  make  use  of  the  fact  that  varying  the  polar  angle  e  relative  to  the 
surface  leads  to  selective  enhancement  of  surface  or  bulk  signals.®  At  a 
given  polar  angle  0,  the  Auger  intensity  from  a  component  characterized  by 
an  atom  number  density  function  p(y)  within  the  overlayer  is  given  by,® 

I  -  aAnD  c  i  n o'*  00 

I(0)  =  fe^‘R^lW^pWexp^y/Asin0^dy'  (1) 

IQ  and  A  are  the  incident  electron  beam  flux  and  cross-sectional  area, 
respecti vely,  0  is  the  Auger  cross  section,  a  is  the  solid  angle  of 
acceptance,  D  is  the  detector  efficiency,  a  is  the  angle  of  incidence  with 
respect  to  the  surface,  R  is  the  electron  reflection  coefficient  for  Auger 
electrons  passing  through  the  surface  potential  barrier,  0"  and  0  are  the 
propagation  angles  inside  and  outside  the  material,  and  A  is  the  electron 
mean  free  path.  The  perpendicular  depth  y  is  measured  relative  to  the 
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surface.  In  order  to  utilize  Eq.  I  to  determine  p(y)  for  components  in  an 
interface,  the  group  of  factors  outside  the  integral  (collectively  referred 
to  hereafter  as  F(e))  must  be  evaluated  separately  for  pure  substrate  and 
overlayer  materials.  This  task  is  readily  accomplished  by  measuring  1(0) 
for  atomically  smooth  surfaces  of  pure  metal  and  semiconductor,  and  solving 
Eq.  1  for  F(©),  recognizing  that  p(y)  is  simply  the  bulk  density  of  the 
material.  Equation  1  can  then  be  evaluated  for  both  metal  and  semiconductor 
in  a  particular  interface  using  F(0)  and  physically  reasonable  choices  for 
p(y).  Theory  can  then  be  fit  to  experiment  to  determine  the  optimal  choice 
for  p(y ) . 

We  have  employed  an  Auger  electron  spectrometer  with  angle-resolving 
capability  described  elsewhere.7  Single  crystal  wafers  of  Si (111)  were 
cleaned  in  situ  by  cycles  of  Ar  ion  bombardment  and  annealing  to  800°C.  The 
resulting  surfaces  showed  a  sharp  reconstructed  7x7  LEED  pattern  and  were 
free  of  contaminants,  as  judged  by  Auger  spectra  collected  at  0=10°  relative 
to  the  surface. 

High  purity  Cu  and  Au  were  separately  evaporated  onto  a  Si  substrate 
using  a  quartz  crystal  oscillator  to  monitor  coverage.  Pressures  never 
exceeded  2x10"*®  Torr  at  any  stage  of  evaporation  or  measurement  (base 
pressure  =  4x10"**  Torr).  Polar  scans  were  constructed  from  spectra  taken 
every  2°  from  0=90°  (normal  emission)  to  0=0°.  LEED  shows  that  at  coverages 
above  ~10A,  the  Cu /Si ( 11 1 )  surface  is  ordered  with  a  (111)  orientation 
rotated  30°  with  respect  to  the  substrate.  Therefore,  polar  scans  in  two 
symmetry  inequivalent  azimuthal  planes  perpendicular  to  the  surface  were 
taken  and  averaged  to  remove  intensity  fluctuations  due  to  diffraction 
effects.  Since  Au  overlayers  are  disordered,  polar  scans  in  only  a  single 
azimuthal  plane  were  collected. 
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In  Figs.  1  and  2,  we  show  Auger  spectra  as  a  function  of  coverage  for 
the  Cu  and  Au  interfaces.  In  the  bottom  spectrum  in  each  figure,  sufficient 
metal  has  been  deposited  to  bury  essentially  all  the  Si,  even  though  at  150A 
a  small  amount  of  Si  diffuses  to  the  surface  of  the  Au  overlayer.  Polar 
profiles  taken  at  these  coverages  are  used  to  determine  F(e)  for  pure  Cu  and 
Au,  although  a  slight  correction  is  necessary  for  the  latter  to  account  for 
the  small  amount  of  Si  still  present  on  the  surface.  From  similar  spectra 
taken  at  normal  emission,  we  have  constructed  attenuation  curves  in  the  form 
of  tn[ISi  (d)/Is^  ( o ) 3  vs.  d,.  where  d  is  the  overlayer  thickness  in  A.  A 
sharp  interface  would  render  this  function  linear,  with  slope  equal  to  -1/A. 
In  both  cases,  the  plots  exhibit  positive  deviations  from  linearity, 
indicating  Si  outdiffusion.  The  initial  slope  provides  a  good  estimate  of 
A.  We  obtain  2.8±.5A  and  5.0±.4A  for  Cu  and  Au,  respectively.  These  values 
are  used  in  the  evaluation  of  Eq.  1.  Both  metals  induce  either  a  splitting 
or  a  broadening  of  the  Si (L2  3VV)  line  above  a  coverage  of  2A,  indicating 
disruption  of  the  Si  lattice  and  modification  of  the  Si  valence  states  as 
interaction  with  the  metal  and  outdiffusion  proceed.®  Qualitatively,  the 
persistence  of  the  Si (Lg  3VV)  line  at  coverages  as  high  as  20A  indicates 
significant  outdiffusion. 

For  both  interface  systems,  the  optimal  choice  of  Si  density  function 
is  given  by 

p0  exp[-a(d-y)]  for  y<d  (2a) 

y  P0  for  y>d  (2b) 

where  p0  is  the  bulk  density  for  Si (.05018  atoms/A®),  d  is  the  overlayer 
thickness,  and  a  is  a  free  parameter  chosen  to  maximize  agreement  between 
theory  and  experiment.  Table  1  shows  optimal  values  of  a  and  the  predicted 
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Si  density  on  the  surface  for  coverages  ranging  from  2  to  20A.  Also,  Fig.  3 
shows  the  experimental  and  theoretical  polar  profiles  based  on  the 
parameters  in  Table  1.  The  Cu-Si  intermixed  phase  is  ordered  and  the  Si 
profiles  at  high  coverages  show  strong  diffraction  modulation,  in  spite  of 
averaging  over  two  azimuthal  planes.  The  disorder  of  the  Au-Si  intermixed 
phase  is  reflected  by  the  lack  of  structure  in  Si  profiles. 

It  is  useful  to  define  a  characteristic  diffusion  length  (CDL)  as  the 
diffusion  depth  of  either  Si  in  the  overlayer,  or  metal  in  the  substrate,  by 
which  the  density  has  dropped  to  1/e  of  its  value  at  the  interface  boundary. 
As  shown  in  Table  I,  the  CDL  increases  gradually  from  18%  of  the  overlayer 
thickness  at  a  Au  coverage  of  2A  to  39%  at  20A,  whereas  in  Cu  the  CDL 
abruptly  increases  from  20%  of  the  overlayer  thickness  at  2A  to  ~35%  at  5, 

10  and  20A. 

For  both  metals,  the  best  choice  of  density  function  is  given 

ps+al-y  for  ancl  a^  coverages  (3a) 

p(y)  ={  (ps+ajdJexpC-agCy-d)]  for  y>d  and  d<10A  (3b) 

not  measurable  for  y>d  and  d=20A  (3c) 

Here,  ps  is  the  metal  density  at  the  surface  and  is  allowed  to  vary  freely. 
When  necessary,  a^  and  a£  were  constrained  to  be  <0  and  >0,  respectively. 

In  this  way,  these  functions  represent  a  metal  density  which  decreases 
linearly  from  surface  to  interface  and  decays  exponentially  in  the 
substrate.  The  maximum  probing  depth  is  ~3a  (~11A  for  both  metal  CVV  Auger 
electrons).  Using  a  technique  described  elsewhere,6  values  of  3.6  and  3.4A 
were  determined  for  the  Cu^^VV)  and  Au(Ng>7VV)  electrons,  respecti vely. 
The  resulting  values  of  p$,  a^  and  ag  are  shown  in  Table  1.  Interestingly, 
ai  converged  to  zero  for  5  and  10A  Cu  coverages  and  would  have  been  positive 
unless  constrained.  Similarly  for  2,  4  and  8A  Au  coverages,  83  converged  to 


zero  and  would  have  been  negative  otherwise.  Despite  the  constraints,  the 
fits  are  quite  good  and,  at  higher  coverages,  the  metal  density  at  the 
surface  is  calculated  to  approach  that  of  pure  metal  for  both  systems. 

Figure  4  shows  the  experimental  and  theoretical  polar  profiles  for  both 
metal  overlayer  systems.  As  in  the  Si  polar  profiles  for  the  well-ordered 
Cu /Si (111)  system,  the  Cu  polar  profiles  show  considerable  diffraction  fine 
structure  at  higher  coverages,  despite  azimuthal  averaging. 

The  results  are  best  visualized  by  plotting  the  calculated  densities  as 
a  function  of  distance  from  the  interface  boundary,  as  shown  in  Figs.  5  and 
6.  The  picture  which  emerges  is  that  Au  reacts  more  vigorously,  initiates 
more  intermixing  with  longer  metal  CDL  values,  and  exhibits  a  gradual  onset 
of  reaction.  Cu,  on  the  other  hand,  appears  to  trigger  Si  outdif fusion  at 
coverages  in  excess  of  2A,  as  witnessed  by  an  order  of  magnitude  increase  in 
surface  Si  density  in  going  from  2A  to  5A  thicknesses,  but  does  not 
penetrate  as  deeply;  CDL  values  are  limited  to  1  to  2A.  These  results  are 
consistent  with  the  extensive  Si  lattice  disruption  brought  about  by  Au  but 
not  by  Cu,  leading  to  a  disordered  overlayer  for  the  former.  Moreover,  the 
persistence  of  surface  Si  even  at  Au  thicknesses  of  150A  is  in  stark 
contrast  with  the  complete  disappearance  of  the  Cu(M2  3VV)  signal  by  a 
coverage  of  60A. 

Through  the  use  of  polar-angle  resolved  Auger  electron  emission  and 
appropriate  modeling,  we  have  extracted  quantitative  information  about 
composition  and  diffusion  phenomena  in  the  Au/Si ( 111 )  and  Cu /Si (111) 
interface  systems.  Such  quantitative  information  is  extremely  important  in 
modeling  quantitative  interface  behavior  and  nicely  complements 
photoemission  spectroscopy,  which  provides  insight  into  the  chemical 
character  of  the  intermixed  phase. 
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TABLE  1 

Model  calculations  of  composition  and  diffusion  phenomena 

based  on  experimental  polar  profiles  1 

Cu/Si (111) 


d  ( A ) 

a 

Si  surface 
density*. 

( atoms /lr) 

Si  CDL 

(A) 

al 

a2 

Cu  surface 

density** 

(atoms/AJ) 

Cu  CDL 
(A) 

2 

2.51 

3.3xl0"4 

-2.9xl0-2 

.667 

l.OxlO"1 

1.5 

5 

2.5xl0“3 

1.7 

0 

.452 

6.9xl0"2 

2.2 

.282 

3.5 

0 

.923 

8.2xl0“2 

1.1 

.143 

2.9xl0~3 

6.7 

3.8xl0"3 

- 

7.9xl0“2 

- 

• 

Au/Si  (111-) 

d(A) 

a 

Si  surface 
density*. 

( atoms/A J) 

Si  COL 
(A) 

al 

a2 

Au  surface 

density* 

(atoms/AJ) 

Au  CDL 
(A)* 

2 

2.78 

1.9xl0-4 

.36 

-5.8xl0-3 

0 

2.9xl0“2 

00 

4 

1.14 

5.3xl0“4 

.88 

-2.7xl0’3 

0 

3.6xl0‘2 

Op 

8 

.345 

3.2xl0“3 

2.9 

-1.9xl0-3 

0 

4.4xl0“2 

OP 

.128 

3.9xl0"3 

7.8 

-2.8xl0"3 

• 

5.7xlO“2 

•m 

*Bulk  Si  density  *  5.0xl0-2  atoms/A3 

**Bulk  Cu  density  =  8.5xl0“2  atoms/A3 
+Bulk  Au  density  *  5.9xl0-2  atoms/A3 


^''Infinite"  means  to  the  limit  of  depth  sensitivity  (3a  or  ~10A) 


9- 


FIGURE  CAPTIONS 

Fig.  1  Auger  spectra  for  the  Cu/Si (111 )-7x7  interface  as  a  function  of 

coverage.  The  incident  electron  beam  energy  and  current  were  4  keV 
and  0.1  to  0.5  uA,  respectively.  Spectra  were  collected  at  a  polar 
angle  of  48°  relative  to  the  surface,  and  0=0°  lies  in  the  (110) 
azimuthal  plane  perpendicular  to  the  surface.  1A  of  Cu  equals  1.1  ml 
on  the  Si (111 )  surface. 

Fig.  2  Auger  spectra  for  the  Au/Si (111 )-7x7  interface  and  a  3  keV  primary 
electron  energy.  1A  of  Au  equals  0.75  ml  on  the  Si (111)  surface. 

Fig.  3  Experimental  and  theoretical  polar  intensity  profiles  for  the  Si (L2  3VV) 
line  in  the  two  interface  systems-  The  top  panels  show  the  smoothed 
clean-surface  polar  profiles  used  in  determining  F(o). 

Fig.  4  Experimental  and  theoretical  polar  intensity  profiles  for  the  two  metal 
CVV  lines.  The  intensities  in  the  top  right-hand  panel  have  been  in¬ 
creased  proportionately  to  stimulate  a  completely  Si -free  Au  surface. 

Fig.  5  Cu  and  Si  atom  number  densities  as  a  function  of  distance  from  the 
interface  boundary,  as  determined  by  fitting  intensities  calculated 
using  Eq.  1  to  those  from  experiment.  Dashed  lines  indicate  extra¬ 
polations  beyond  the  maximum  probing  depth  of  3a. 

Fig.  6  Au  and  Si  atom  number  densities  as  a  function  of  distance  from  the 
interface  boundary,  as  determined  by  fitting  intensities  calculated 
using  Eq.  1  to  those  from  experiment.  Dashed  lines  indicate  extra¬ 
polations  beyond  the  maximum  probing  depth  of  3a. 
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Modeling  A  Heterogeneous  Metal/Semiconductor  Interface:  Ce  on  Si (111) 

M.  Grioni,  J.  Joyce*,  M.  del  Giudice,  D.G.  O'Neill*  and  J.H.  Weaver 

Department  of  Chemical  Engineering  and  Materials  Science 
University  of  Minnesota,  Minneapolis,  MN  55455 

Abstract 

High  resolution  synchrotron  radiation  photoemission  studies  of  Ce 
deposited  onto  cleaved  Si ( 111 ) -2x1  reveal  heterogeneous  growth  which 
involves  clustering,  Ce/Si  reaction  to  form  silicide  patches,  lateral 
silicide  growth,  and  finally  Ce  overlayer  formation  with  surface  segregated 
Si,  Core  level  lineshape  analysis  reveals  three  distinct. Si  local  bonding 
configurations.  The  relative  interface  concentration  of  each  Si  specie  has 
been  determined  as  a  function  of  overlayer  thickness,  and  a  model  for  this 
interface  is  presented  and  discussed. 
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Interface  interactions  are  often  modeled  by  uniform  layers  of  reacted 

or  unreacted  material  on  stable  substrates.  Although  the  assumption  of  a 

homogenous  interface  simplifies  modeling,*  it  is  not  correct  for  many 
2  16 

interfaces.  Indeed,  interface  reactivity,  kinetics,  and  chemical 

trapping  are  difficult  to  treat  within  the  homogeneous  layer  model. 

Further,  treatments  of  Schottky  barrier  formation  based  on  uni  form 

overlayers  fall  short  of  adequately  describing  the  behavior  of  the  barrier 
2  3 

at  low  coverage.  Unfortunately,  while  the  importance  of  heterogeneity 
has  long  been  recognized,  few  studies  have  been  able  to  resolve  lateral 
features  at  an  interface. 

In  this  paper,  we  discuss  an  interface  having  a  richness  of  phenomena 
related  to  heterogeneous  growth:  Ce/Si (111).  As  we  will  show,  results 
obtained  with  high  resolution  core  level  photoemission  make  it  possible  to 
model  even  complicated  systems  by  distinguishing  the  different  electronic 
configurations  of  the  constituents—  Such  core  level  studies  are 
complemented  by  valence  band  results  which  emphasize  changes  in  the  states 
near  the  Fermi  level  Ep. 

In  a  previous  paper, we  discussed  the  submonolayer  behavior  of  the 
Ce/Si (111)  interface,  demonstrating  cluster  formation  by  combining  core 
level  and  valence  band  photoemission,  LEED,  and  angle-resolved  Auger 
spectroscopy.  Here,  we  emphasize  the  development  of  the  interface 
subsequent  to  the  ripening  of  the  clusters  (reaction),  describing  the 
varying  environments  which  evolve  with  coverage  at  room  temperature. 

The  photoemission  experiments  were  conducted  at  the  Wisconsin 
Synchrotron  Radiation  Center  using  the  Tantalus  storage  ring  and  the 
grasshopper  and  toroidal  grating  monochromators.  The  overall  resolution  of 
the  measurements  ranged  from  200  meV  at  hv=40  eV  to  500  meV  at  hv=135  eV. 


Ce  was  deposited  onto  cleaved  Si ( 111 ) -2x1  surfaces  in  an  experimental  system 
described  in  detail  elsewhere.  During  evaporation,  the  pressure  rose  from 
operating  pressures  of  3x10"^  to  ~lxl0"^  Torr.  Evaporation  rates  of 
~lA/min  and  the  use  of  a  shutter  permitted  0.04  ml  incremental  depositions 
without  difficulty.  In  this  paper,  coverages  will  be  defined  in  monolayers 
where  1  ml =2.6  A  =  7.8xl0iH  atoms -cm  ,  the  surface  density  of  Si. 

In  Fig.  1  we  show  valence  band  energy  distribution  curves  (EDCs)  for  Ce 
coverages  0-16  ml.  Although  the  EDCs  have  been  drawn  with  approximately  the 
same  height,  comparison  of  absolute  count  rates  shows  that  the  Ce-derived 
features  very  quickly  dominate  the  Si  emission,  indicative  of  the  higher 
photoionization  cross  section  of  Ce  5d4f  states  relative  to  Si  sp  states  at 

IQ 

hv=40  eV.  The  difference  curve  at  the  bottom  of  Fig.  1  gives  the  result 
of  subtraction  of  normalized  EDCs  for  clean  Si  and  0.4  ml  Ce/Si.  For  this 
difference  curve,  the  well-defined  doublet  1.2  and  3.2  eV  below  Ep  reveals 
emission  from  Ce-induced  states  ancLis  indicative  of  weak.ly  interacting  Ce 
clusters,  as  proven  by  valence  band,  LEED  and  angle-resol ved  Auger  results 
(Ref.  17).  With  increasing  coverage,  the  spectral  features  change  as  peak 
"A"  near  Ep  is  lost  and  the  valence  band  becomes  relatively  structureless. 

At  that  point,  the  interface  photoemission  spectra  can  no  longer  be 
described  as  the  superposition  of  unreacted  Si  and  Ce  clusters.  By  0=2.4  ml 
the  interface  product  is  metallic  with  a  distinct  Fermi  level  cutoff.  At  4 
ml  coverage,  a  sharp  feature  emerges  at  Ep  and  grows  relative  to  the  deeper 
structure  and  the  valence  band  narrows,  ultimately  converging  to  Ce.  Figure 
1  therefore  shows  that  valence  band  emission  can  be  readily  identified  at 
low  coverage  (clusters)  and  high  coverage  (Ce  metal  overlayer),  but 
continuous  evolution  at  intermediate  coverage  suggests  a  broad,  intermixed 
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interface  or  a  heterogeneous  interface.  Core  level  studies  make  it  possible 
to  rule  out  the  former. 

The  results  on  the  left  of  Fig.  2  show  the  Si  2p  core  emission  measured 

with  high  surface  sensitivity  (hv=135  eV, escape  depth  ~4A^).  The 

bottom-most  EDC  for  clean  Si (lll)-2xl  reveals  the  characteristic  2p  doublet 

on 

broadened  by  surface-shifted  components.  At  submonolayer  coverage,  it 
shifts  rigidly  to  greater  binding  energy  (band  bending  without  chemical 
reaction).17  At  9=0. 6ml  a  reacted  component  appears  and  grows  relative  to 
the  original  component,  as  indicated  by  the  tic  marks  shifted  0.67  eV.  For 
0>3ml ,  a  third  Si  doublet  is  observed  at  1.2  eV  lower  binding  energy  than 
bulk  Si.  This  component  persists  to  very  high  coverage  and  the  Si  2p 
lineshape  is  significantly  sharper  than  for  the  clean  surface,  suggesting  a 
more  atomic  Si  bonding  configuration. 

On  the  right  of  Fig.  2,  we  show  the  deconvolution  of  the  experimentally 
observed  Si  2p  emission  into  three .spin-orbit -split  pairs,  with  branching 
ratios  and  spin-orbit  splittings  equal  to  those  of  Ref.  20.  For  each,  the 
experimental  full  width  at  half  maximum  (FWHM)  was  derived  from  results  at 
9=0.2  ml  where  the  surface  shifted  component  was  small  and  the  linewidth  was 
a  measure  of  our  resolution.  Comparison  shows  this  FWHM  to  be  ~0.1  eV 
larger  than  for  the  clean  surface,  analogous  to  what  has  been  observed  by 
Ludeke  et  al  in  detailed  examinations  of  the  nonreactive  Ag/GaAs(110) 
interface.  Such  broadening  is  physically  reasonable  for  the  room 
temperature  reacted  phase  because  of  the  likely  existence  of 
not-quite-identical  bonding  configurations  for  Si,  i.e.  to  disorder  within 
the  dominant  local  bonding  configuration.  The  significant  conclusion  from 
Fig.  2  is  that  the  experimental  lineshapes  can  be  fit  very  well  for  all 
coverages  by  assuming  only  three  unique  Si  configurations  (three  doublets). 
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The  attenuation  of  the  total  Si  2p  emission  is  a  measure  of  the  rate  at 
which  Si  is  masked  by  the  overlayer.  However,  for  a  reacted  interface  with 
several  components,  it  is  much  more  informative  to  consider  the  attenuation 
of  each  component  separately.  In  Fig.  3  we  display  the  attenuation  results 
for  the  three  Si  components,  Oj  =  ln[I(2pi (o)/I(2p(initial )].  As  shown, 
substrate  Si-1  attenuation  is  rapid  throughout  the  coverage  range.  The 
reacted  Si-2  component  grows  with  coverage  from  0.6  ml  to  about  2.5  ml,  but 
then  diminishes  sharply.  The  Si -3  component  appears  slightly  above  3  ml  and 
remains  nearly  constant  with  coverage,  amounting  to  about  8%  of  the  original 
Si  intensity  at  0=24  ml  (not  shown). 

An  interface  with  behavior  demonstrated  by  Figs.  1-3  is  complex,  but 
the  results  presented  here  allow  formulation  of  a  model  which  describes  its 
heterogeneous  growth.  Although  clusters  form  below  0.6  ml,  they  interact 
weakly  with  the  substrate  and  the  Si  2p  binding  energy  does  not  change.  At 
0.6  ml,  conversion  from  cluster  growth  to  reactive  intermixing  occurs. 
Nonetheless,  upon  reaction  the  disrupted  surface  area  is  still  only  a 
portion  of  the  total  surface,  having  resulted  from  intermixing  beneath  and 
adjacent  to  the  clusters.  The  emergence  of  the  Si-2  signal  reveals  Si  in 
these  patches  and  comparison  with  unreacted  Si  shows  the  chemical 
inequivalence  of  the  two  configurations. 

At  0.6  ml  coverage,  where  the  reacted  Si  component  is  first  observed, 
the  featureless  valence  band  spectra  give  no  evidence  for  a  metallic  Ce 
layer.  However,  both  the  valence  bands  and  the  core  levels  do  indicate  that 
a  Ce  layer  forms  shortly  thereafter.  In  the  valence  bands,  this  can  be  seen 
by  the  reappearance  nkar  Ep  of  peak  "A"  -  which  ultimately  becomes 
recognizable  as  Ce  metal  at  high  coverage.  In  the  core  levels,  the  fact 
that  Si -2  reaches  a  maximum  near  2.5  ml  and  is  then  rapidly  attenuated 


indicates  a  covering-up  of  the  reacted  specie  for  0>2.5  ml.  Nonetheless, 
the  rate  of  attenuation  of  this  reacted  Si -2  specie  is  nearly  exponential 
with  a  characteristic  length  of  -6A,  i.  e.  an  escape  depth  larger  than 
expected  for  a  uniform  overlayer  of  an  efficient  scatterer  like  Ce. 

Likewise,  the  attenuation  of  the  unreacted  Si-1  specie  suggests  the  presence 
of  uncovered  Si  substrate  at  relatively  high  coverage. 

The  results  for  the  intermediate  coverage  regime  can  be  reconciled  by 
modeling  the  interface  as  a  patchwork  of  reacted  Si  on  a  matrix  of  unreacted 
substrate.  Ce  deposited  onto  the  clean  matrix  will  have  the  high  mobility 
of  Ce  observed  at  lowest  coverage  and  can  diffuse  to  the  edge  of  the  reacted 
patches.  The  perimeter  of  these  patches  provides  active  sites  for  continued 
reaction  due  to  the  disruption  of  the  substrate  covalent  bonds.  As  a 
result,  Ce  silicide  growth  {i.  e.  conversion  of  Si-1  to  Si-2)  involves  a 
portion  of  Ce  adatoms  in  the  coverage  range  0.6-3  ml  and  the  patches  grow 
laterally.  This  continues  until  mort  of  the  surface  is  reacted  at  about  3 
ml  nominal  coverage. 

In  the  intermediate  coverage  range,  Ce  deposition  onto  the  reacted 
patches  must  be  treated  differently  from  that  onto  the  unreacted  Si  surface, 
as  indicated  by  the  behavior  of  the  Si-2'  specie.  The  results  of  Figs.  1-3 
indicate  that,  although  the  patches  grow  laterally,  they  apparently  do  not 
grow  vertically.  Instead,  the  reacted  layer  is  rapidly  covered  up  by  Ce,  as 
shown  for  example  by  the  onset  of  the  Ce-like  valence  band  features  and  the 
decay  of  the  Si -2  reacted  specie.  The  slow  Si -2  attenuation  can  be  seen  to 
result  from  the  balance  between  Si -2  attenuation  by  the  overlayer  and  Si -2 
creation  by  the  laterally-growing  patches. 


Deconvolution  of  the  core  emission  demonstrates  that  the  third  Si 
specie  appears  at  nominal  coverages  of  3  ml  and  all  three  Si  components  are 
present  in  about  equal  strength  in  the  surface  region  for  for  0=5  ml  (~8%  of 
the  original  signal).  At  this  coverage,  the  valence  band  results  show  a 
metallic  Ce  overlayer,  and  the  Si -2  attenuation  reveals  a  covering-up  of  the 
reacted  phase.  Indeed,  the  appearance  of  Si -3  is  coincident  with  the 
emergence  of  peak  A  in  the  valence  bands.  This  strongly  suggests  that  Si -3 
is  associated  with  the  Ce  film..  Its  persistence  and  nearly  constant 
concentration  suggests  that  Si -3  is  surface  segregated  on  the  Ce  film.  For 
coverages  greater  than  5  ml,  this  floating  Si  is  the  most  prevalent  specie. 
By  8-10  ml  there  is  no  effectively  no  unreacted  Si-1  or  reacted  Si -2  near 
the  surface,  consistent  with  the  completion  of  silicide  formation  and  Ce 
covering-up.  Indeed,  by  this  coverage  the  valence  bands  are  converging  on 
bulk  Ce.  The  almost  complete  loss  of  Si-1  and  Si-2  signals  by  8  ml 
indicates  that  the  Ce-Si  reaction  i_s  confined  to  a  thin,  well-defined  layer. 
In  the  inset  of  Fig.  2  we  schematically  show  our  model  of  the  evolving  Ce/Si 
interface,  indicating  the  lateral  growth  of  the  silicide  phase,  its 
covering-up  and  the  appearance  of  Si  on  the  Ce  overlayer. 

Comparison  of  the  present  results  with  other  rare-earth/Si  interfaces 
show  similarities  and  intriguing  differences.  All  the  systems  studied  so 
far  (Ce/Si, 17  Sm/Si,22  Eu/Si,22  Yb/Si2^)  exhibit  a  critical  coverage  for  the 
reaction  and  large  chemical  shifts  indicative  of  substantial  charge  transfer 
in  the  reacted  phase.  On  the  other  hand,  while  the  Ce/Si  interface  is 
characterized  by  a  low  value  of  the  Schottky  barrier,  as  expected  on  the 
basis  of  recent  results  for  technological  samples,12  no  such  lowering  of  the 
barrier  was  found  for  the  Sm/Si  system.  This  difference  is  not  understood 
at  present  and  a  systematic  study  of  Schottky  barrier  formation  involving 
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rare-earth  metals  and  Si  is  needed.  Moreover,  the  fact  that  the  chemistry 
of  rare-earth/Si  systems  is  dominated  by  a  critical  temperature*^  suggests 
that  interesting  new  effects  in  their  evolution  could  be  expected  as  a 
function  of  temperature.  Work  in  this  direction,  particularly  concerning 
the  possible  existence  of  ordered  phases  at  high  temperature,  is  currently 
in  progress. 

Finally,  a  significant  point  of  this  paper  has  been  that  high 
resolution  studies  of  reactive  interfaces  make  it  possible  to  distinguish 
the  chemical  environments  of  the  constituents,  i.  e.  that  specific  bonding 
configurations  form  which  can  be  identified  even  for  intermixed  interfaces 
where  long  range  order  is  absent  and  the  scale  of  the  interface  is  limited 
to  a  few  monolayers.  Hence,  we  can  conclude  that  strong  local  bonding 
determines  the  character  of  the  evolving  interface  species.  Comparison  to 
bulk  compounds  will  show  how  the  interface  phase  fits  into  the  hierarchy  of 
the  bulk  phase  diagram.  For  Ce/Si,_we  find  no  evidence  that  more  than  a 
single  silicide  phase  exists,  and  whatever  stoichiometry  gradients  might 
exist  between  the  reacted  patches  and  the  unreacted  Si  are  too  small  for 
identification.  By  studying  the  variation  of  each  component  with  nominal 
metal  overlayer  coverage,  even  complex  interfaces  can  be  modeled.  The 
implications  of  these  conclusions  are  major  because  they  indicate  the 
limitations  of  modeling  of  interface  properties  based  on  homogeneous 
overlayers  and  they  show  the  importance  of  microscopic  fingerprinting  of 
interface  phases. 
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Figure  Captions 

FIG.  1.  Energy  distribution  curves  for  Ce/Si  showing  valence  band 
evolution.  The  bottom-most  curve  represents  the  normalized  difference  between 
clean  Si (111  )-2xl  and  the  interface  with  0.4  ml  coverage.  At  low  coverage  the 
spectra  are  superpositions  of  Ce  and  Si,  at  high  coverage  they  converge  to 
bulk  Ce,  but  in  the  intermediate  range,  they  reveal  a  heterogeneous  interface. 

FIG.  2.  EDCs  on  the  le*t  show  the  evolution  of  the  Si  2p  core  emission. 
On  the  right,  lineshape  deconvolutions  are  shown  for  representati ve  coverages. 
At  any  coverage,  the  results  can  be  fit  with  components  representing  clean  Si 
(Si-1),  a  reacted  phase  (Si-2),  and  surface  segregated  Si  (Si-3).  Our  model 


n 


for  the  interface  is  shown  pictorally  in  the  inset. 

FIG.  3.  Attenuation  of  the  total  Si  2p  emission  (solid  line  at  top)  and 
the  attenuation  of  each  of  its  three  components.  Although  the  substrate 
attenuation  is  rapid  and  monotonic,  the  reacted  component  grows  to  a  maximum 
near  2.5  ml  and  is  then  attenuated.  The  surface-segregated  component  appears 
only  after  the  reacted  component  is  covered-up  by  Ce  and  is  nearly  constant  in 
magnitude  to  high  coverage. 
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ABSTRACT. 

Si(L2,3VV)  and  Cu(M2j3VV)  Auger  intensities  from  the  Cu/Si ( 111 )-7x7 
interface  have  been  measured  in  a  polar-angle  resolved  fashion  for  various 
Cu  thicknesses.  A  simple  theory  of  polar-angle  resolved  Auger  emission  has 
been  developed  and  used  to  extract  atom  number  densities  for  both  Si  and  Cu 
as  a  function  of  depth.  It  is  found  that  significant  Si  outdiffusion  is 
triggered  at  300K  by  Cu  coverage  in  excess  of  2A.  Cu  diffuses  into  the 
substrate  with  a  characteristic  penetration  depth  of  1  to  2A  at  all 
coverages  and,  above  a  coverage  of  2A,  Si  diffuses  into  the  overlayer  with  a 
characteristic  depth  of  -35%  of  the  overlayer  thickness.  The  Cu  density 
appears  to  be  constant  throughout  the  intermixed  region. 


PACS:  64.75.+g 
66.30. Ny 
73.40.-c 
79.60.Gs 
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JL  INTRODUCTION 

Condensation  of  metal  atoms  on  single-crystal  semiconductor  surfaces 
leads  to  a  wide  variety  of  interesting  physical  and  chemical  phenomena.1 
Ordered  layer-by-layer  growth,  cluster  formation, ^  interdiffusion, and 
chemical  reaction®  have  been  observed  for  different  systems.  Successful 
observation  of  these  phenomena  requires  a  wide  repertoire  of  experimental 
probes;  no  one  technique  can  provide  all  the  necessary  information  to 
unambiguously  model  the  microscopic  development  of  the  Interface. 

In  this  paper  we  address  the  issue  of  interdiffusion  at  a 
metal -semi conductor  interface.  This  research  was  undertaken  because 
knowledge  of  the  atomic  densities  of  the  substrate  and  adsorbed  material  as 
a  function  of  depth  in  the  Interface  region  is  critical  for  characterizing 
interface  properties.  Host  attempts  to  obtain  such  Information  have 
utilized  argon  Ion  sputter  profiling.6"11  As  discussed  by  some  of  these 
authors,  °  ion-induced  surface  damage,  preferential  sputter  rates, 
implantation  of  sputtered  species,  depth  calibration  problems,  and  poor 
depth  resolution  render  the  quantitative  reliability  of  this  technique 
dubious,  particularly  when  applied  to  ultra-thin  overlayers  (a  few 
angstroms).  In  order  to  overcome  these  difficulties,  attempts  have  been 
made  to  obtain  a  better  understanding  of  the  sputtering  process®*11  and  to 
quantitatively  account  for  the  phenomenon  of  sputter  broadening.® 

A  nondestructive  alternative  to  sputter  profiling  which  eliminates  many 
of  the  above  problems  utilizes  the  polar  angle  dependence  of  Auger  and 
photoelectrons  emitted  from  both  the  substrate  and  overlayer.  As  originally 
demonstrated  by  Fadley  and  Bergstrom, varying  the  collection  angle 
relative  to  the  plane  of  the  surface  leads  to  selective  enhancement  of  bulk 
and  surface  signals.  Most  such  studies  to  date  have  utilized  core-level 
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photoelectrons  excited  by  either  conventional  x-ray  sources'1,  or  synchrotron 
radiation. ^  One  disadvantage  of  polar-angle  resolved  photoemission  is  that 
of  data  acquisition  time.  In  order  to  obtain  a  complete  polar  intensity 
profile,  a  large  number  of  spectra  must  be  taken  over  the  collection  angle 
range  0  to  90°.  A  minimum  of  nine  or  ten  spectra  would  allow  a  rough  polar 
profile  to  be  constructed,  but  20  to  40  spectra  would  be  desirable.  Active 
metal  overlayers  will  contaminate  in  the  several-hour  time  span  required  to 
acquire  such  detailed  data  by  means  of  photoemission. 

A  closely  related  yet  much  faster  approach  to  the  problem  involves 
utilizing  the  polar  angle  dependence  of  low-energy  core-valence-valence 
(CVV)  Auger  electron  Intensity.  An  excellent  spectrum  can  be  obtained  in 
the  pulse-counting  N (E )  mode  in  just  a  few  minutes.  Moreover,  low  incident 
beam  currents  (a  fraction  of  a  uA)  can  be  used  to  minimize  surface  damage 
and  one  still  obtains  very  high  count  rates.  In  this  article  we  describe 
the  first  use  of  polar-angle  resolved  Auger  emission  to  quantitatively  study 
Interdiffusion  phenomena  at  a  metal -semiconductor  interface.  The 
Cu/SI (111)— 7x7  system  was  chosen  because  Interdiffusion  of  Si  into  the 
overlayer  has  been  observed  and  qualitatively  studied  with  other  surface 
analysis  techniques, including  XPS,  UPS,  AES,  and  LEED. 

This  article  is  organized  as  follows:  Section  II  describes  the  theory 
of  polar-angle  resolved  Auger  emission  as  developed  for  this  particular 
application.  Section  III  provides  experimental  details.  Section  IV  presents 
experimental  and  theoretical  results,  Section  V  consists  of  a  discussion  of 
the  work,  and  conclusions  are  drawn  in  Section  VI. 


II.  THEORY  OF  POLAR-ANGLE  RESOLVED  AUGER  EMISSION  FROM  INTERMIXED  OVERLAYERS 

After  an  Auger  electron  has  been  ejected  from  an  atom  in  a  solid,  it 
can  be:  (1)  inelastical ly  scattered,  leading  to  signal  attenuation,  (2) 
elastically  scattered,  leading  to  diffraction  effects  if  the  sample  is  a 
single  crystal,  and  (3)  refracted  at  the  surface-vacuum  interface.  In  this 
model  we  treat  the  over layer  as  a  uniform  slab  with  a  constant  attenuation 
coefficient  and  ignore  any  single-crystal  effects. 

Figure  1  shows  the  relevant  geometry  of  our  sample-analyzer 
arrangement.  Auger  electrons  originating  a  perpendicular  distance  y  from 
the  surface  propagate  with  internal  angle  e'  to  refract  at  the  surface  and 
travel  toward  the  analyzer  aperture  with  external  angle  0.  The  principal 
factor  which  governs  intensity  variations  with  depth  is  inelastic 
attenuation.  Thus,  the  intensity  at  a  given  depth  y  will  be  proportional  to 
exp(-y/Asine' ),  where  A  is  the  electron  mean  free  path.  .  In  addition, 
surface  refraction  will  bring  about- changes  in  intensity  due  to  internal 
reflection  and  what  is  effectively  an  increase  in  Auger  electron  flux  caused 
by  a  reduction  in  the  cross  sectional  area  of  the  Auger  electron  beam  upon 
refraction.  These  effects  add  factors  of  ( 1-R )  and  sine '/sine, 
respectively,  where  R  is  a  reflection  coefficient  for  Auger  electrons  within 
the  sample  incident  on  the  surface.^  Additional  factors  governing  the 
observed  intensity  are  incident  electron  beam  flux  (IQ),  Auger  cross  section 
(a),  sample  area  Irradiated  (A/s1na),  solid  angle  of  acceptance  (fl), 
detector  efficiency  (D),  and  emitting  atom  number  density  at  depth  y  (p(y)). 
Combining  these  factors,  the  differential  intensity  brought  about  by  a 
volume  element  a  perpendicular  distance  y  from  the  surface  is  given  by 
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dl(e,y)  *  [-^7^-][(l-R)(^^-)]p(y)exp(-y/Asine,)dy.  (1) 

The  propagation  angle  within  the  material  ( ©' )  is  related  to  that 
outside  by  the  conservation  of  k,  rule  which  yields, 

9'  ■  cos~*[(-!L-— cose].  (2) 

fck 


Ek  is  the  electron  kinetic  energy  inside  the  material  and  VQ  is  the  inner 
.potential.  These  refraction  effects  are  very  important  at  the  kinetic 
energies  of  the  CVV  Auger  transitions  used  here  (-100  eV),  particularly  at 
low  emission  angles.  As  shown  by  Fadley,15  the  electron  reflection 
coefficient  can  be  estimated  by  treating  the  interaction  of  the  electron 
with  the  surface  as  analogous  to  a  particle  penetrating  a  square  well 
potential  of  height  VQt  the  Inner  potential  of  the  material.  From  this 
model  we  get 

R  •  [(l-(l-V0/Eks1n20')1/2)/(l+(l-V0/Elts1n29')1/2)]2.  (3) 

With  the  exception  of  the  cross  section,  the  factors  within  the  first 
set  of  brackets  in  Eq.  1  depend  only  on  the  spectrometer.  The  factors  in 
the  second  set  of  brackets,  which  account  for  the  change  in  Auger  electron 
flux  brought  about  by  surface  refraction,  can  be  combined  with  the 
instrument  factors  and  cross  section  to  yield  a  function  F(e),  defined  as: 


(A) 
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This  function  accounts  for  everything  except  the  emitting  atom  number 
density  in  the  material  and  electron  attenuation,  and  should  be 
approximately  constant  as  coverage  proceeds,  provided  the  inner  potential 
and  cross  section  do  not  change  substantially.  Thus,  if  F(e)  is  evaluated 
for  each  component  in  the  interface,  it  could  be  used  to  deduce  atomic 
densities  in  the  intermixed  phase.  For  single-crystals  of  Cu  and  Si,  the 
inner  potentials  are  14  eV  and  17  eV,  respectively. Therefore,  we  take 
15  eV  as  an  approximate  value  for  the  Cu-Si  interface  at  all  stages  of 
development.  In  practice,  we  find  that  the  calculated  intensities  are  not 
particularly  sensitive  to  the  choice  of  VQ  over  the  range  14  to  17  eV. 
Changes  in  CVV  Auger  cross  section  for  either  Cu  or  Si  would  result  from 
valence  band  modification  upon  intermixing  and  are  revealed  by  lineshape 
changes.  No  such  changes  have  been  observed  for  the  Cu^^VV)  line.  The  - 
S1(L2,3W)  line  does  split  and  broaden  at  higher  Cu  coverages  and  such, 
effects  are  not  incorporated  into  the  model. 

For  pure  samples  of  either  substrate  or  overlayer  material,  p(y)  is 
simply  equal  to  p,  the  bulk  atom  number  density  for  each  material.  In  this 
case,  Eq.  (1)  Is  easily  integrated  over  y  and  F (e)  can  be  directly 
evaluated.  The  result  is 


F(e) 


1(9) 

pAsine’ 


(5) 


Once  determined  for  pure  substrate  (s)  and  overlayer  (o)  materials,  F$(e)  - 
and  F0O)  will  be  Inserted  into  integrated  forms  of  Eq.  1  appropriate  to  the 


interface. 
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1(e),  and  hence  F(e),  as  used  In  Eq.  5  represent  an  average  over  the  8° 
polar  angle  resolution  of  our  analyzer.  For  later  use  it  is  desirable  to 
express  F(e)  in  units  of  per  degree  of  polar  angle.  In  order  to  make  this 
adjustment,  it  is  assumed  that  for  all  polar  angles  except  88°-90°  and  0-2°, 
1(e)  increases  linearly  with  e  over  all  8°  polar  angle  intervals.  This 
approximation  is  reasonable  in  light  of  the  observed  behavior  of  1(e)  for 
both  substrate  and  overlayer  materials  (to  be  discussed  in  more  detail  jn~- 
Section  IV  A).  It  then  suffices  simply  to  divide  1(e)  by  8°  for  each  value 
of  e  prior  to  evaluating  Eq.  5.  The  fact  that  1(e)  is  observed  to  go 
through  a  maximum  at  0*90°  and  go  to  zero  for  0<O°  is  used  to  perform  a 
suitable  adjustment  to  obtain  an  1(0)  per  degree  of  polar  angle  for  0*88-90° 
and  0-2°. 

In  order  to  extract  concentration  information  from  the  interface  polar 
profiles,  a  particular  form  of  p(y)  for  both  substrate  and  overlayer 
material  Is  assumed.  Equation  (1)  can  then  be  integrated  to  yield 
expressions  for  I$(0)  and  I0( ©) .  For  the  substrate  material, 

Is(0)  *  Fs(0)/ps(y)exp(-y/Assin0,s)dy  (6) 

and  similarly  for  the  overlayer.  The  function  p(y)  contains  one  to  three 
undetermined  parameters  which  are  chosen  to  generate  the  best  fit  with 
experiment..  The  fitting  algorithm  used  is  a  combination  of  gradient  search 
and  linearization  of  the  fitting  function  methods  developed  by  Marquart.18 
Further  details  on  the  exact  choice  of  p(y)  are  given  in  Section  IV  A. 

Finally,  comparison  of  theoretical  results  with  experimental  polar 


profiles  requires  that  the  former,  which  come  out  per  degree  of  polar  angle, 
be  Integrated  over  the  finite  aperture  size  of  8°.  This  task  is  readily 
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accomplished  using  Simpson's  and  trapezoidal  rules  in  the  calculation  of 
1(9). 

III.  EXPERIMENTAL  DETAILS 

The  system  used  to  perform  angle-resolved  Auger  spectroscopy  is 
described  in  detail  elsewhere.19  The  single-pass  CMA  with  angle-resolving 
capability  is  also  used  as  a  LEED  I-V  detector,  emplo  ig  a  modification 
first  described  by  McDavid  and  Fain.^®  Single-crystal  wafers  of  (111) 
oriented  p-type  Si  were  cleaned  by  cycles  of  Ar  ion  sputtering  and  heating 
to  800°C.  This  procedure  produced  a  clean  7x7  surface;  no  impurities  were 
'present  as  judged  by  Auger  spectra  taken  at  a  surface-enhancing  polar  angle 
of  10°. 

Evaporation  of  high-purity  Cu  was  monitored  by  a  quartz  crystal 
oscillator.  Pressures  during  evaporation  never  exceeded  2xl0"10  Torr  and 
during  polar  scans  were  typically  ?-9xl0-11  Torr.  Complete  polar  scans  from 
normal  emission  to  grazing  emission  were  obtained  by  taking  spectra  every 
2°.  Replicate  measurements  at  fixed  angles  indicated  a  standard  deviation 
of  3-4%  of  the  value  of  the  mean.  The  total  time  required  to  complete  a 
fulT  polar  scan  was  about  two  hours  and  grazing  emission  Auger  spectra  taken 
at  the  end  of  each  run  showed  that  the  surface  was  free  of  contaminants. 
Polar  profiles  were  obtained  in  the  (110)  azimuthal  plane  perpendicular  to 
the  surface  and  in  a  plane  rotated  30°  from  (110).  These  were  then  averaged 
to  remove  as  much  diffraction  modulation  as  possible,  since  our  theoretical 
model  ignores  single-crystal  effects.  For  future  reference,  we  define  the 
zero  in  i  as  lying  in  the  (110)  plane  perpendicular  to  the  surface. 
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4  KeV  electrons  were  used  to  excite  the  Auger  transitions  and  incident 
beam  currents  were  typically  Q.l-0.5uA.  There  was  no  evidence  of 
electron-induced  surface  damage  at  any  point,  as  judged  by  reproducible 
Auger  intensities  and  LEED  patterns  over  several  hours  of  beam  Irradiation. 
Peak  intensities  were  determined  by  smoothing  the  spectra,  substracting  a 
linear  background,  and  Integrating.  Spectra  were  then  normalized  to  one 
another  by  dividing  by  the  average  number  of  background  counts  on  the  flat 
high-energy  side  of  the  Si (L2  3VV)  peak,  thus  removing  Intensity  variations 
brought  about  by  drifts  in  the  Incident  beam  current. 

IV.  RESULTS 

Figure  2  demonstrates  how  the  Auger  spectrum  develops  as  a  function  of 
Cu  coverage.  Coverages  are  expressed  In  Angstroms,  and  one  Angstrom  of  Cu 
equals  1.1ML  on  the  Si ( 111 )  surface.  From  0-2A  there  1s.no  significant 
change  In  the  Si (L£  3VV)  lineshape.-  However,  above  5A  a  triplet  develops, 
indicating  substantial  modification  of  the  valence  states  of  Si  as  the 
normal  sp3  hybrid  bonds  are  broken  and  Si  diffuses  into  the  Cu  layer. 

-  Similar  Hneshape  changes  have  been  observed  as  a  function  of  coverage  for 
Pt/SI (lll),^-  Pd/S1(Ill),^  and  Cu/Si (111 ) .^  Using  spectra  taken  at  normal 
emission,  one  can  construct  a  SiCLg^W)  attenuation  curve  to  determine  the 
extent  to  which  interdiffusion  occurs. 

Figure  3  shows  a  plot  of  An ( I (d)/I (0))  vs  d  where  d,  1(0),  and  1(d)  are 
the  overlayer  thickness,  clean  surface  S1(L£  3VV)  intensity,  and  SKL2  3VV) 
intensity  with  overlayer  thickness  d  present,  respectively.  Each  point  is 
an  average  of  four  spectra  taken  at  normal  emission  and  the  standard 
deviation  of  each  data  set  is  less  than  4%  of  the  mean  for  that  set.  Linear 
behavior  with  slope  -1/A  is  expected  if  the  interface  is  sharp,  and  the 
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observed  positive  deviation  from  linearity  indicates  outdif fusion  of  Si. 

The  initial  slope  provides  a  good  estimate  of  the  mean  free  path  for  the 
Si^^W)  electron  moving  through  the  overlayer.  In  this  case,  we  get 
2.8±.5A.  In  the  evaluation  of  Eq.  5,  we  use  a  value  of  3A.  Also  shown  for 
reference  is  a  growth  curve  for  the  Cu(M2  3VV)  line. 

Corresponding  LEED  data  indicate  relaxation  of  Si(lll)-7x7  to  a  rapidly 
weakening  Si(lll)-lxl  pattern  for  2  and  5  A  coverages?-  fol lowed  by  the 
appearance  at  10A  of  a  Cu(lll)-lxl  pattern  rotated  30°  with  respect  to  the 
Si  surface  mesh.  These  results  are  in  good  agreement  with  those  obtained  by 
Ringeisen  et  al .3 

In  Figs.  4  and  5  we  present  polar  angle  intensity  profiles  for  the 
Si(^2,3vV)  and  Cuff^^W)  Auger  intensities  in  two  symmetry -inequivalent 
azimuthal  planes.  The  top  profiles  in  each  panel  correspond  to  pure  Si 
(Fig.  4)  or  Cu  (Fig.  5)  and  contain  the  information  from  which  F$(e)  and 
FQ(e)  are  evaluated  (Eq.  5).  At  a  coverage  of  60A,  no  Si (*-2,3vv )  signal  was 
present  and  the  LEED  I-V  curves  were  those  expected  for  a  Cu(lll) 
single-crystal  surface  film.  Thus  we  take  this  surface  as  our  reference 
surface  for  pure  Cu(lll). 

The  diffraction  modulation  observed  in  each  case  is  a  sensitive  probe 
of  the  local  structure  of  the  emitting  atom  and  can  be  used  as  a  structural 
fingerprint  at  other  coverages.  At  2A  coverage,  the  diffraction  modulation 
observed  for  the  clean  Si (111) -7x7  suface  is  weakened  but  still  present  and 
the  Cu  profile  Is  essentially  featureless.  Coupled  with  the  weak  Si (111)1x1 
LEED  pattern  and  the  attenuation  curve  in  Fig.  3,  these  data  suggest  that 
the  Cu  layer  grows  in  a  weakly  ordered  layer-by-layer  fashion,  up  to  2A. 

This  mode  of  growth  diminishes  and  gives  rise  to  reaction  and  outdif fusion 
at  5A,  as  evidenced  by  the  Si  lineshape  change  and  deviations  from  linearity 
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of  the  $i(l-2  3W)  attenuation  curve  (Fig.  3).  By  10A,  the  Cu  polar  profiles 
show  clear  similarities  to  those  of  the  60A  film  (pure  Cu ( 111 ) ) .  However,  a 
significant  amount  of  Si  remains  dispersed  throughout  the  overlayer  and  the 
increased  Si(L£  3VV)  splitting  indicates  a  more  substantial  chemical 
interaction  with  the  Cu.  There  are  no  substantive  changes  between  10A  and 
20A,  except  that  the  overlayer  is  only  slightly  richer  in  Cu  but  poorer  in 
Si  by  a  factor  of  2.  At  60A,  the  evolution  of  the  overlayer  to  a  pure 
Cu(lll)  film  is  complete  and  no  Si  is  present  on  the  surface. 

In  the  top  panels  of  Figs.  6  and  7,  we  present  average  polar  profiles 
for  clean  Si (111 )— 7x7  and  a  60A  Cu  overlayer  (essentially  a  pure  Cu(lll) 
surface).  Each  profile  is  an  average  over-two  symmetry  inequivalent 
azimuthal  planes  (d=o°  and  d=30°).  As  can  be  seen,  such  averaging  is  not 
sufficient  to  remove  all  the  diffraction  features.  Therefore,  for  the 
purpose  of  evaluating  F$(e)  and  FQ(9),  we  effectively  smooth  out  the 
remaining  diffraction  modulation.  Then,  after  adjusting  the  intensities  for 
8°  angular  resolution  as  discussed  in  Section  II,  these  curves  are  used  to 
evaluate  F(e)  for  Si(l-2,3VV)  and  Cu^^VV)  Auger  emission  through  the  use 
of  Eq.  5. 

Also  shown  in  Figs.  6  and  7  are  the  experimental  polar  profiles 
averaged  over  azimuthal  planes  at  «J*o°  and  d=30°  compared  to  those  predicted 
by  Eq.  6  for  the  best  choice  of  density  function.  For  both  S i (L£  3VV)  and 
Cy(M2,3vV)  emission  calculations,  physically  plausible  atom  number  density 
functions  of  various  forms  with  one,  two,  or  three  undetermined  parameters 
were  tried.  For  each  function,  the  fitting  routine  determined  the  choice  of 
parameters  which  yielded  a  minimum  value  of  x^.  Turning  first  to  the 
distribution  of  Si  in  the  interface,  a  linear  function  and  an  exponential 
function  were  modeled.  The  linear  function  was  of  the  form 
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(7a) 

(7b) 


p(>) 


^p0  -  a(d-y)  for  y<d 

p0  for  y>d 


and  the  exponential  function  was  of  the  form 

.  p0exp[-a(d-y )]  for  y«d  (8a) 

p  (y  )  a  I  u 

p0  for  y>d.  (8b) 

Here  p0  is  the  bulk  atom  number  density  of  Si (5.018xl0“3  atoms/A3),  d  is  the 

overlayer  thickness,  a  is  the  undetermined  parameter,  and  y  is  the 

perpendicular  distance  measured  from  the  surface.  It  is  found  that  the 

.exponential  function  .gives  a  much  better  fit  to  the  experimental  data  than 

does  the  linear  function;  x  for  the  exponential  function  was  1.5  times 

lower  at  2A  coverage  and  3  times  lower  at  20 A  coverage  than  that  for  the 

— - 

linear  function.  For  both  functions,  a  mean  free  path  of  3A  was  used,  as 
discussed  earlier.  The  optimal  values  of  a  and  the  predicted  Si  atom  number 
densities  on  the  surface  are  given  in  Table  I. 

Clearly,  the  calculated  surface'Si  density  is  an  order  of  magnitude 
lower  at  2A  than  that  for  the  5,'T07_and  20A  overlayers,  indicating  a 
substantial  increase  in  outdiffusion  above  2A.  Moreover,  for  all  coverages 
above  2A,  the  characterfstfc  penetration  depth  (distance  into  the  overlayer 
'by  which  the  Si  density  has  decayed  to  1/e  of  its  bulk  value)  is 
approximately  35%  of  the  overlayer  thickness,  whereas  for  2A  it  is  only 
about  20%  of  the  overlayer  thickness.  These  results  are  supported  by  the 
observed  change  in  Si (L2 ^ 3VV )  lineshape  above  2A. 

For  the  Cu  distribution  calculation,  we  have  used  a  linear  function 
which  terminates  di scontinuously  at  the  Cu-Si  interface  and  a  linear 
function  with  an  exponential  tail  into  the  substrate.  In  analytical  form, 
the  functions  are 


and 


(9a) 

(9b) 


p(y)  -  { 


0S+aly  for  y<cl 
0  for  y>d 


P(y)  .  {°S  +  *1*  for  y<d 

(ps+a1d)exp[-a2(y-d)]  for  y>d. 


(10a) 

(10b) 


Here,  ps  is  the  Cu  number  density  at  the  surface  (a  free  parameter),  d  and  y 
are  as  defined  previously,  and  aj  and  a2  are  undetermined  parameters.  The 
-Cu(M2  3VV)  mean  free  path  was  determined  by  performing  a  two-parameter  fit 
on  the  20A  coverage  and  employing  a  linear  density  function  of  the  form 
p(y)a. 078+ajy,  where  a^  and  A  were  allowed  to  vary.  The  value  .078  was 
chosen  because  values  of  .078  to  .079  routinely  result  from  two-parameter 
fits  using  the  form  p(y)=ps+aiy  with  a  fixed  value  of  A  taken  from  the 
universal  mean  free  path  curve.  This  calculation  generates  a  value  of 

3.6A  for  A,  which  is  used  in  all  subsequent  calculations^  In  general,  the 

2  ~ 

X  values  are  slightly  lower  for  the  linear  functions  (9a&b)  than  for  the 
linear-exponential  functions  (10a4b)  for  all  coverages.  However,  from  the 
point  of  view  of  using  physically  reasonable  results,  we  judge  the 
linear-exponential  functions  (lOa&b)  to  be  most  appropriate  for  the  2,  5, 
and  10A  coverages  and  the  linear  functions  (9a&b)  to  be  best  for  the  20A 
coverage. 

Qualitative  considerations  of  the  2  and  5A  polar  profiles  suggest  that 
if  Cu  grows  epitaxially  with  no  indiffusion,  the  Cu(M2  3VV)  intensity  should 
be  larger  at  shallow  exit  angles,  by  virtue  of  surface  enhancement. 

However,  the  intensity  is  essentially  constant  from  -30°  to  90°  for  both 
coverages.  This  effect  could  be  brought  about  by  a  substantial  increase  in 
Cu  density  at  the  interface  boundary  relative  to  the  surface,  or  by  Cu 
diffusion  into  the  substrate.  Either  phenomenon  will  bolster  the 
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intensities  at  large  polar  angles.  These  two  situations  are  predicted  to 
occur  with  the  use  of  Eqs.  9aAb  and  lOaAb,  respectively.  If  the  system  is 
constrained  to  produce  the  observed  polar  profile  without  the  aid  of  metal 
indiffusion,  a  large  increase  in  Cu  density  at  the  interface  relative  to  the 
surface  must  result.  Equations  9a&b  predict  increases  from  .0542  and  .0528 
(at  the  surface)  to  .1342  and  .1161  atoms/A3  (at  the  interface)  for  2A  and 

.  5A  coverages,  respectively.  However,  the  weak  lxl  LEED  patterns  at  2  and  5A 
suggest  that  within  the  first  few  monolayers  the  Cu  density  should  be  close 
to  that  of  Si  (.0502  atoms/A  ).  Moreover,  given  that  the  Si  number  density 
increases  from  surface  to  interface  boundary,  it  is  most  likely  that  the  Cu 
number  density  decreases  or  remains  constant.  Therefore,  we  judge  Eqs.  9a&b 
to  be  physically  untenable  representations  of  the  Cu  density  function.  In 
contrast,  Eqs.  IGa&b  predict  for  2A  a  sharp  decrease  in  Cu  density  from  .101 
at  the  surface  to  .0424  atoms/A3  at  the  interface.  Using  Eqs.  lOaAb  for  the  - 
5A  coverage,  it  was  necessary  to  constrain  a^  to  be  negative  and  less  than 
or  equal  to  ps/d  in  magnitude.  Without  this  constraint,  a^  tended  to  be 
positive,  as  it  was  when  Eqs.  9a£b  were  used.  With  the  constraint,  however, 
aj  converages  to  zero,  resulting  in  a  constant  Cu  density  of  .0687  atoms/A3 
throughout  the  interfacial  region.  For  both  coverages,  the  extent  of  Cu 
indiffusion  is  predicted  to  be  about  the  same;  characteristic  penetration 
depths  of  1  to  2A  result  from  the  fits. 

In  principle,  Eqs.  lOa&b  can  be  used  for  the  10A  coverage  as  well. 

-Since  95%  of  a  given  Auger  signal  originates  within  a  depth  of  3a  ( 11 A  for 
the  Cu(M£  3VV)  electron),  we  can  probe  the  interface  boundary  when  d  is  ~1QA 
or  less.  As  in  the  5A  case,  aj  must  be  bracketed  between  zero  and  -ps/d, 
lest  it  converge  to  a  positive  number.  The  result  is  a  constant  Cu  density 
of  .0817  atoms/AJ  within  the  overlayer  and  a  characteristic  penetration 
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depth  of  1A  in  the  substrate.  The  use  of  Eqs.  9a£b  without  constraints  on 
a^  results  in  an  unreasonably  large  and  positive  value  of  a^  (2.067x10“^). 
Since  we  have  no  sensitivity  to  the  interface  boundary  when  d  is  20A,  we  use 
Eqs.  9a£b  to  calculate  Cu  densities  in  the  top  10  to  11  A  of  the  overlayer. 
Here,  we  do  not  constrain  aj  and  a  slight  Cu  density  increase  is  predicted 
from  surface  to  interface  boundary.  Moreover,  the  Cu  density  converges  to  a 
value  quite  close  to  the  bulk  Cu  value  at  the  surface,  in  good  agreement 
with  the  LEED  results  which  showed  the  ordered  Cu ( 111 )  pattern. 

The  overall  results  are  best  suinnarized  by  plotting  the  calculated  Si 
and  Cu  atom  number  densities  as  a  function  of  distance  from  the  interface 
boundary,  as  shown  in  Fig.  8.  One  striking  observation  is  that  the  Si 
surface  density  is  an  order  of  magnitude  smaller  and  the  Cu  surface  density 
is  20  to  30%  larger  at  2A  than  at  higher  coverages.  This  result  indicates 
that  Cu  builds  up  on  the  surface  prior  to  the  strong  Si-Cu  interaction,  the 
surface  disruption,  and  the  onset  of  Si  outdiffusion.  This  interface 
therefore  exhibits  a  reaction  triggered  by  Cu  coverage  in  excess  of  2A  but 
not  at  lower  coverage.  A  second  important  observation  is  that  for  the  10 
and  20A  coverages,  the  total  atom  number  density  within  the  first  few 
Angstroms  of  the  interfacial  region  exceeds  that  of  either  pure  Si  or  pure 
Cu  but  rapidly  approaches  the  bulk  Cu  value  with  increasing  distance  from 
the  interface  boundary.  This  result  suggests  that  the  interface  is  an 
intermixed  Cu-Si  phase,  with  Cu  being  the  predominant  element  and  Si  being 
dissolved  in  the  matrix  to  a  decreasing  extent  as  the  surface  is  approached. 


Interdiffusion  at  metal -semiconductor  interfaces  has  been  qualitatively 
studied  recently  by  both  angle-resolved  and  angle-integrated 
photoemissionl*2-6*24*26.  outdiffusion  of  Si  is  relatively  easy  to 

observe  by  constructing  attenuation  curves  similar  to  those  in  Fig.  3, 
diffusion  of  the  metal  into  the  substrate  has  been  observed  by  comparing  the 
rate  of  growth  of  the  metal  photoemission  intensity  with  coverage  to  that 
expected  for  a  sharp  interface24  and  by  performing  marker  experiments.26  In 
the  latter,  a  thin  layer  of  an  inert  metal  is  deposited  onto  the  substrate 
prior  to  evaporation  of  the  metal  of  interest.  Assuming  the  marker  layer 
does  not  interfere  with  interdiffusion,  the  behavior  of  the  substrate  to 
marker  intensity  ratio  with  coverage  Indicates  qualitati vely  whether  or  not 
Indiffusion  of  the  metal  Is  occurring.  The  advantages  of  the  method 
described  here  Is  that  no  marker  layer  (which  may  indeed  perturb  the 
interdiffusion  process26)  is  needed  and  the  quantitative  extent  of 
Indiffusion  can  be  obtained. 

Other  quantative  methods  Involving  polar  angle-resolved  x-ray 
photoemission  have  been  proposed  recently.  Specifically,  Pijolet  and 
Hoi  linger  have  attempted  to  extract  density  information  directly  from  the 
polar  profiles  via  a  simplex  method  which  minimizes  the  difference  between 
theoretical  and  experimental  intensity  r  tios.2^  However,  this  approach  has 
proved  to  possess  some  serious  difficulties.  First,  the  extracted 
concentrations  are  extremely  sensitive  to  the  details  of  the  polar  profile. 
Intensity  ratio  changes  on  the  order  of  0.001%  make  substantial  changes  in 
the  extracted  concentrations.  In  order  to  remedy  the  situation,  they  found 
that  it  was  necessary  to  introduce  two  physically  reasonable  contraints  on 
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the  concentrations.  These  consisted  of:  (1)  bracketing  the  concentration  of 
each  specie  between  zero  and  the  bulk  concentration  for  the  pure  material 
and  (2)  assuming  a  particular  shape  or  trend  in  the  density  profile. 

Clearly,  once  such  constraints  are  introduced,  the  method  becomes  quite 
analogous  to  what  we  present  here.  However,  this  "direct  extraction"  method 
requires  a  number  of  involved  numerical  techniques  which  necessitate 
considerable  computer  time  on  a  large  mainframe  computer.  The  method  we 
introduce  is  relatively  simple  and  can  be  run  on  any  of  the  powerful 
microcomputers  currently  available. 

B.  The  Chemical  Nature  of  the  Interface  Region 

The  substantial  lineshape  changes  in  the  SifL^^VV)  peak  clearly 
demonstrate  that,  upon  interdiffusion,  a  strong  electronic  interaction 
occurs  between  SI  and  Cu.  As  the  normal  sp3  hybrid  bonds  are  broken  in  the  - 
Si  lattice  and  outdif fusion  occurs,  the  valence  states  are  modified  and  new 
channels  for  Auger  de-excitation  are  opened,  thus  modifying  the  SUl^^W) 
spectrum.^  Although  the  details  of  chemical  environment  are  very  difficult 
to  extract  from  Auger  lineshape  analysis,  the  results  presented  here  enable 
us  to  comment  on  the  stoichiometry  of  the  intermixed  phase.  Above  a 
coverage  of  2A,  the  Cu  density  in  the  overlayer  is  constant.  However,  the 
Si  density  drops  very  rapidly  from  interface  to  surface  for  all  coverages. 
These  results  suggest  that  a  copper  silicide  of  fixed  stoichiometry  does  not 
form  in  the  Intermixed  region,  as  proposed  earlier.4  Rather,  a  Cu  lattice 
which  is  rich  in  Si  only  near  the  interface  appears  to  be  a  better 
description  (see  Fig.  8).  This  conclusion  is  consistent  with  valence  band 
photoemission  results  which  show  a  dominant  Cu-derived  3d  feature  that 
shifts  from  3.5  eV  below  the  Fermi  level  at  submonolayer  coverages  to 


-2.8  eV  at  45  monolayers.  Indicative  of  a  transition  from  isolated  Cu  atoms 
to  a  pure  cu  metallic  phased  Cu  diffuses  into  the  substrate  with  a 
characteristic  penetration  depth  of  1-2A  and,  in  all  likelihood,  occupies 
defect  and  interstitial  sites.  This  description  of  the  interface  is 
consistent  with  the  bulk  Cu-Si  phase  diagram  which  shows  -2%  solubility  of 
Si  in  Cu  at  room  temperature  and  no  compound  formation  at  room  temperature 
and  Si  percent  by  weight  less  than  -8%.^8*^9  The  lack  of  solubility  of  Cu 
in  Si,  however,  indicates  that  the  slight  indiffusion  of  Cu  observed  here  is 
the  result  of  microscopic  localized  phenomena. 

Although  similar  agreement  between  interface  interdiffusion  results  and 
bulk  phase  information  has  been  noted  by  Brillson  et  al_.^8  for 
Al/Si ( 111 )-7x7,  agreement  for  Au/Si (111 )-7x7  is  lacking.  The  former  is  a 
weakly  interacting  chemisorption  system,  consistent  with  phase  data.  The 
latter,  however,  strongly  interacts  with  considerable  Intermixing  at  room 
temperature,  despite  the  inertness  predicted  by  the  phase  diagram.  Thus,  it 
appears  that  bulk  thermodynamic  data  are  not  of  universal  utility  in 
predicting  interface  phenomena  for  noble  metals  adsorbed  on  Si  single 
crystals. 


VI.  CONCLUSIONS 

Angle-resolved  Auger  electron  emission  has  been  shown  to  be  a  powerful 
technique  for  the  study  of  interdiffusion  at  a  metal -semiconductor 
interface.  Detailed  antitative  information  on  the  composition  of  the 
Interfacial  region  can  be  extracted  from  the  polar-angle  intensity  profiles 
by  means  of  a  relatively  simple  theoretical  model.  It  has  been  shown  that 
tntermixing  for  Cu/Si (111 )-7x7  leads  to  an  interface  which  is  constant  in  Cu 
density  throughout  but  which  exhibits  a  Si  density  that  decreases 


exponentially  from  interface  boundary  to  surface.  Furthermore,  it  has  been 
demonstrated  that  the  bulk  phase  diagram  for  the  noble  metal -semiconductor 
system  is  useful  in  predicting  interdiffusion  properties  of  the  interface, 
although  not  universally  so. 
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TABLE  1 


Model  predictions  for  surface  s’ ’■’con  number  densities 


PSf (y)  *  pQexp[-a(d-y )] 


d  (A) 

a 

surface  atom  number 

2 

2.51 

3.34x10 

5 

.602 

2.47x10 

10 

.282 

3.00x10 

20 

.143 

2.85x10 

-3 


Table  2 


Model  predictions  for  Cu  atom  number  densities 


P2+a1y  for  y<d.at  all  coverages 
PCu(y)  *  {  0  Tor  y>d  and  d»2,5  and  10 A 

(ps+aid)expC-a2(y-d)]  for  y>d  and  d«20A 


d(A) 

p„  (surface  atom  number  density  in  atoms/A**) 

a^xlO2) 

a? 

2 

■  -  .101 

-2.93 

.667 

5 

.0687 

0 

.452 

10 

.0817 

0 

.923 

20 

.0793 

.0127 

- 
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FIGURE  CAPTIONS 

Fig.  1  Orientation  of  the  sample  with  respect  to  the  incident  electron  beam  and 
analyzer  aperture. 

Fig.  2  Si  (L2. 3VV )  and  Cu(M2  3W)  Auger  spectra  as  a  function  of  Cu  coverage  in 
A.  1A  of  Cu  equals  1.1  monolayers  on  Si (111). 

'Fig.  3  Attenuation-growth  curve  for  the  unanriealed  Si-Cu  interface,  d  is  the 
overlayer  thickness. 

Fig.  4  Si (L2  3VV)  polar  intensity  profile  taken  in  two  symmetry -inequivalent 
azimuthal  planes  perpendicular  to  the  surface.  No  annealing. 

Fig.  5  Cu (M2  3VV)  polar  intensity  profile  taken  in  two  symmetry -inequi valent 
azimuthal  planes  perpendicular  to  the  surface.  No  annealing. 

Fig.  6a  Clean  surface  SHl^^VV)  polar  profile  averaged  over  two  azimuthal 
planes,  raw  data  and  smoothed  to  remove  diffraction-induced  fine 
structure. 

Fig.  6b  Two-azimuth  averaged  experimental  and  calculated  Si (L2  3VV)  polar 
profiles  using  Eqs.  8a4b  as  a  Si  atom  number  density  function. 

Fig.  7a  Cu(M2  3VV)  polar  profiles  for  a  60A  Cu  film  (pure  Cu(lll))  averaged 

over  two  azimuthal  planes,  raw  data  and  smoothed  to  remove  diffraction- 
induced  fine  structure. 

Fig.  7b  Two-azimuth  averaged  experimental  and  calculated  Cu(M2  3W)  polar 
profiles  using  Eqs.  lOa&b  for  2,  5  and  10  A  coverages  and  Eqs.  9a 
and  9b  for  the  20  A  coverage. 

Fig.  8  Calculated  Cu  and  Si  atom  number  densities  vs.  distance  from  the  inter¬ 
face  boundary.  A  dashed  line  signifies  an  extrapolation,  since  our 
maximum  depth  sensitivity  is  3a  (10  to  11  A). 
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